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A  wide  range  of  practical  problems  can  be  simulated  using  these  shapes. 
For  example,  flat  plates  can  be  used  to  model  the  superstructure 
of  a  ship,  the  body  of  a  truck,  or  the  wings  and  stores  of  an  aircraft. 
The  finite  elliptic  cylinder  can  be  used  to  model  a  mast  or  smoke 
stack  of  a  ship,  or  the  fuselage  and  engines  of  an  aircraft. 

This  document  is  designed  to  give  an  overall  view  of  the  operation 
of  the  computer  code,  to  instruct  a  user  in  how  to  model  structures, 
and  to  show  the  validity  of  the  code  by  comparing  various  computed 
results  against  measured  data  whenever  available.  It,  also,  describes 
in  detail  the  input  and  output  data  for  the  code.  This  information 
should  be  sufficient  for  most  user's  to  learn  how  to  effectively 
operate  the  code. 
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CHAPTER  I 
INTRODUCTION 


The  Numerical  Electromagnetic  Code  -  Basic  Scattering  Code  is 
a  user-oriented  computer  code  for  the  analysis  of  the  far  field  pat¬ 
terns  of  antennas  in  the  presence  of  perfectly  conducting  metal  struc¬ 
tures  at  UHF  and  above.  Complicated  structures  can  be  simulated 
by  arbitrarily  oriented  flat  plates,  an  infinite  ground  plane,  and 
a  finite  elliptic  cylinder.  This  type  of  analysis  has  been  used 
very  successfully  in  the  past  to  model  aircraft  shapes  [1,2,3 ].  The 
present  solution  has  been  extended  to  include  a  wide  range  of  problems. 
For  example,  flat  plates  can  be  used  to  model  the  superstructure 
of  a  ship,  the  body  of  a  truck,  or  the  wings  and  stores  of  an  aircraft. 
The  finite  elliptic  cylinder  can  be  used  to  model  a  mast  or  smoke 
stack  of  a  ship,  or  the  fuselage  and  engines  of  an  aircraft. 

The  analysis  is  based  on  uniform  asymptotic  techniques  formu¬ 
lated  in  terms  of  the  Geometrical  Theory  of  Diffraction  (GTD)  [4,5,6 1. 

The  GTD  approach  is  ideal  for  a  general  high  frequency  study  of  antennas 
in  a  complex  environment  in  that  only  the  most  basic  structural  fea¬ 
tures  of  an  otherwise  very  complicated  structure  need  to  be  modeled. 

This  is  because  ray  optical  techniques  are  used  to  determine  com¬ 
ponents  of  the  field  incident  on  and  diffracted  by  the  various  sttuc- 
tures.  Components  of  the  diffracted  fields  are  found  using  the  GTD 
solutions  in  terms  of  the  individual  rays  which  are  summed  with  the 
geometrical  optics  terms  in  the  far  field.  The  rays  from  a  given 
scatterer  tend  to  interact  with  other  structures  causing  various 
higher-order  terms.  In  this  way  one  can  trace  out  the  various  possible 
combinations  of  rays  that  interact  between  scatterers  and  determine 
and  include  only  the  dominant  terms.  Thus,  one  need  only  be  con¬ 
cerned  with  the  important  scattering  components  and  neglect  all  other 
higher-order  terms.  This  method  leads  to  accurate  and  efficient 
computer  codes  that  can  be  systematically  written  and  tested.  Com¬ 
plex  problems  can  be  built  up  from  simpler  problems  in  manageable 
pieces. 
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’In'  limitations  associated  with  the  computer  code  result  from 
the  basic  nature  of  the  analyses.  The  solution  is  derived  using 
the  GTD  which  is  a  high  frequency  approach.  In  terms  of  the  scattering 
from  plate  structures  this  means  that  each  plate  should  have  edges 
at  least  a  wavelength  long.  In  terms  of  the  cylinder  structure  its 
major  and  minor  radii  and  length  should  be  a  wavelength  in  extent. 

In  addition,  each  antenna  element  should  be  at  least  a  wavelength 
from  all  edges  and  the  curved  surface.  In  many  cases,  the  wavelength 
limit  can  be  reduced  to  a  quarter  wavelength  for  engineering  purposes. 

Modeling  small  structures  and  antennas  can  be  better  accomplished 
using  an  integral  equation  solution  such  as  NEC-Moment  Methods[71. 

The  Basic  Scattering  Code  has  been  interfaced  with  the  Moment  Method 
code  so  that  the  capabilities  of  both  methods  can  be  used  to  the 
fullest.  For  example,  the  Moment  Method  code  can  be  used  to  analyze 
the  currents  and  impedance  of  an  antenna.  The  magnitude  and  phase 
of  the  current  weights  can  then  be  used  in  the  Basic  Scattering  Code 
to  predict  the  far  field  patterns  of  the  antennas  in  arbitrary  pattern 
cuts. 

There  are  two  documents  describing  the  NEC-Basic  Scattering 
Code.  The  present  document  in  known  as  Part  I.  Part  II  is  a  Code 
Manual [81  that  describes  the  FORTRAN  coding  in  detail.  The  Code 
Manual,  first,  gives  background  on  practical  aspects  of  the  GTD. 

Several  examples  are  shown  to  illustrate  how  the  various  GTD  fields 
superimpose  to  give  a  total  solution.  Next,  a  particular  GTD  term 
is  discussed  in  more  detail  to  show  the  general  concepts  involved 
throughout  the  code.  An  overview  on  how  the  code  is  organized  is 
discussed  along  with  a  description  of  the  various  coordinate  systems 
involved,  how  a  general  subroutine  is  organized,  and  how  the  various 
subroutines  are  interrelated.  The  Code  Manual  also  contains  for 
each  subroutine:  (1)  a  statement  of  purpose,  (2)  an  illustration 
showing  the  geometry  involved,  13)  a  brief  narrative  on  the  method 


used,  (4)  a  flow  diagram,  (5)  a  dictionary  of  major  variables,  and 
(6)  a  listing  of  the  code.  Finally,  it  defines  the  common  blocks 
and  lists  the  system  library  functions  used  by  the  code.  The  in¬ 
formation  in  the  Code  Manual  will  be  of  primary  interest  to  someone 
attempting  to  modify  the  code.  It  will  also  be  helpful  when  the 
code  is  being  implemented  on  a  computer  system  on  which  the  coding 
may  not  be  compatible. 

This  document  is  designed  to  give  an  overall  view  of  the  operation 
of  the  computer  code,  to  instruct  a  user  in  how  to  use  it  to  model 
structures,  and  to  show  the  validity  of  the  code  by  comparing  various 
computed  results  against  measured  data  whenever  available.  Chapter 
II  describes  an  overall  view  of  the  organization  of  the  code.  A 
detailed  description  of  the  input  command  words  and  their  associated 
input  parameters  is  given  in  Chapter  III.  How  to  apply  the  capabilites 
of  this  input  data  to  a  practical  structure  is  briefly  discussed 
in  Chapter  IV.  This  includes  a  clarification  of  the  subtle  points 
of  interpreting  the  input  data.  The  representation  of  the  output 
is  discussed  in  Chapter  V.  Various  sample  problems  are  presented 
in  Chapter  VI  to  illustrate  the  operation,  versatility,  and  validity 
of  the  code.  Most  users  of  the  code  should  find  that  the  User's 
Manual  is  sufficient  to  learn  how  to  effectively  operate  the  code. 
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f'HAPTFR  II 

PRINCIPLES  OF  OPERATION 


T^e  NEC-Basic  Scatte>  ing  Code  is  designed  to  hp  a  user  oriented 
computer  redo.  T^e  rerpssary  data  to  describe  a  p^oMem  can  he  input 
and  the  ••'esiilting  answers  can  he  obtained  with  a  minimum  amount  of 
knowledge  by  the  user  on  how  the  code  operates.  As  with  most  codes, 
however,  it  is  necessary  to  have  at  least  a  basic  knowledge  of  the 
key  points  in  order  to  be  able  to  intelligently  use  it  and  interpret 
its  results.  This  section  is  designed  to  give  just  a  brief  descripti 
of  the  code  for  this  purpose.  The  Code  Manual,  which  is  the  second 
part  of  the  Basic  Scattering  Code  documentation,  gives  more  in  depth 
information  about  the  FORTRAN  coding.  Thus,  this  information  will 
not  be  repeated  here. 

The  NEC -Basic  Scattprinq  Code  is  constructed  in  a  systematic 
wav,  such  that  the  various  ODPrations  of  the  code  are  set  up  in  modu¬ 
lar  spctions.  The  flow  diaqram  shown  helow  illustrates  t^e  major 
divisions  of  th»  main  program.  The  fi*-st  Da>*t  of  tuP  main  nronram 
is  the  input  section  whore  the  geometry  of  the  problem  is  described. 
The  method  used  to  input  data  into  the  computer  code  is  ha^e^  on 
a  command  word  system.  Details  of  the  available  commands  and  the 
ways  to  use  them  are  given  in  the  next  chapter. 

Once  the  necessary  information  to  descrihp  a  Problem  iS  input 
into  the  code,  the  program  analyzes  the  data  and  puts  it  into  the 
correct  form  so  that  the  electric  fields  can  be  calculated.  This 
includes  normalizing  the  geometry  to  wavelengths,  organizing  t^e 
data  into  the  optimum  coordinate  system  for  computations,  and  de¬ 
fining  th*1  fixed  geometry  bounds  for  a  given  source.  Of  course, 
all  of  these  operations  and  the  ones  to  follow  are  done  ODaque  to 
thp  user. 


FLOW  DIAGRAM  OF  THE  MAIN  PROGRAM 
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(Return  to  read  new  command) 
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Table  1 

List  of  Field  Subroutines 


Plate  Field  Subroutines  K=1 


J=1  INCFLD  - 
J=?  REFPLA  - 
J-T  RPLRPL  - 
,1-A  DIFPLT  - 
,l=r'  RPLDPL  - 

J=6  DPLRPL  - 


direct  field 

field  reflected  from  a  plate 

field  doubly  reflected  by  plates 

field  diffracted  by  a  plate 

field  reflected  hv  a  plate  then  diffracted  by 

a  plate 

field  diffracted  by  a  olate  then  reflected  by 
a  plate 


Cylinder  Field  Subroutines  K=? 

J=1  SCTCYL  -  field  scattered  by  a  cylinder 

J=?  REFCAP  -  field  reflected  by  an  end  cap 

J=3  ENDIF  -  field  diffracted  by  an  end  cap  rim 


Plate-Cylinder  Interaction  Field  Subroutines  K=3 

J=1  RPLSCL  -  field  reflected  by  a  plate  then  scattered  by 
a  cylinder 

field  scattered  by  a  cylinder  then  reflected 
by  a  plate 

field  reflected  by  a  cylinder  then  diffracted 
by  a  plate 

field  diffracted  by  a  plate  then  reflected  by 
a  cylinder. 


J=2  SCLRPL  - 
J=3  RCLOPL  - 
J=A  DPLRCL  - 


The  scattering  code  then  computes  the  electric  fields  for  each 
individual  source  in  succession.  Each  GTD  scattered  field  type  is 
hroken  up  into  a  separate  subroutine.  As  can  be  seen  from  the  flow 
chart,  the  code  is  structured  so  that  all  of  one  type  of  scattered 
field  is  computed  at  one  time  for  the  complete  pattern  cut  so  that 
the  amount  of  core  swapping  is  minimized  thereby  reducinq  overlaying 
and  increasing  efficiency.  This  also  is  an  important  feature  that 
allows  the  code  to  be  used  on  small  computers  that  are  not  large 
enough  to  accept  the  entire  code  at  one  time.  The  code  can  be  broken 
into  smaller  overlay  segments  which  will  individually  fit  in  the 
machine.  The  results  are,  then,  superimposed  in  the  main  program 
as  the  various  segments  are  executed. 
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Hip  field  computation  part  of  the  code  is  divided  into  three 
large  sections  ( the  K=1,?,.T  loop!.  The  first  section  (K=!l  contains 
the  major  scattered  fields  associated  with  the  individual  flat  plates 
and  the  interactions  hPtween  the  different  plates.  These  include 
the  direct  field,  the  sinqly  reflected  fields,  douMv  reflec^er1  fields, 
the  sinqlv  diffracted  fields,  the  ref lected-diffracted  fields,  and 
the  diffracte^-ref lect.ed  fields.  The  diffracted  fields  include  the 
normal  diffracted  fields  as  well  as  slope  diffraction,  a  newlv  de¬ 
veloped  heuristic  corner  diffracted  field  and  slope-corner  diffracted 
field.  The  double  diffracted  fields  are  not  included  at  present, 
but  a  warninq  is  provided  wherever  this  field  component  might  be 
important.  This  is  usually  only  a  small  angular  section  of  space. 

Tin's  field  may  be  included  later  whenever  time  and  effort  permit. 

The  second  section  f K=2 )  contains  the  major-  scattered  fields  associated 
with  the  finite  elliptic  cylinder.  This  includes  the  direct  field, 
if  not  already  computed  in  the  plate  section,  the  reflected  field 
the  transition  field,  the  deep  shadow  fields,  the  reflected  field 
from  the  end  caps,  and  the  diffracted  field  from  the  end  cap  rims. 

The  diffracted  field  from  the  end  cap  rim  is  not  at  present  corrected 
in  the  pseudo  caustic  regions.  This  is  where  three  diffraction  points 
on  the  rim  coalesce  into  one.  This  is  only  important  in  small  angular 
regions  in  space  and  is  not  deemed  appropriate  to  he  included  at 
fhe  present  tim».  An  equivalent  current  mpthod  could  be  used  for 
this  small  reoion  but  it  is  rather  time  consuming  to  use  for  the 
^nefits  derived  from  it  for  such  a  qeneral  code.  The  third  serf'pn 
f K = 3 T  contains  the  major  scattered  fields  associated  with  the  inter¬ 
actions  between  the  plates  and  cylinder.  This  includes,  at  present, 
the  fields  reflected  from  the  plates  then  reflected  or  diffracted 
from  the  cylinder,  the  fields  reflected  from  the  cylinder  then  re¬ 
flected  from  the  plates,  the  fields  reflected  from  the  cylinder  then 
diffracted  from  the  plates,  and  the  fields  diffracted  f-'orn  the  plates 
then  reflected  from  the  cylinder.  These  terms  have  been  found  to 
he  sufficient  for  engineering  purposes  when  analyzing  most  structures. 


The  subroutines  for  each  of  the  scattered  field  components  are 
all  structured  in  the  same  basic  way.  First,  the  ray  path  is  traced 
backward  from  the  chosen  observation  direction  to  a  particular  scat- 
terer  and  subsequently  to  the  source  using  either  the  lews  of  reflection 
or  diffraction.  Each  ray  path,  assuming  one  is  possible,  is  then 
checked  to  see  if  it  is  shadowed  by  any  structure  along  the  complete 
ray  path.  If  it  is  shadowed  the  field  is  not  computed  and  the  code 
proceeds  to  the  next  scatterer  or  observation  direction.  If  the 
path  is  not  interrupted  the  scattered  field  is  computed  using  the 
appropriate  GTD  solutions.  The  fields  are  then  superimposed  in  the 
main  program.  This  shadowing  process  is  often  speeded  up  by  making 
various  decisions  based  on  bounds  associated  with  the  geometry  of 
the  structure.  This  type  of  knowledge  is  used  wherever  possible. 

The  shadowing  of  rays  is  a  very  important  part  of  the  scattering 
code.  It  is  obvious  that  this  approach  should  lead  to  various  dis¬ 
continuities  in  the  resulting  pattern.  However,  the  GTD  diffraction 
coefficients  are  designed  to  smooth  out  the  discontinuities  in  the 
fields  such  that  a  continuous  field  is  obtained.  When  a  scattered 
field  is  not  included  in  the  result,  therefore,  the  lack  of  its  pres¬ 
ence  is  apparent.  This  can  be  used  to  advantage  in  analyzing  com¬ 
plicated  problems.  Obviously  in  a  complex  problem  not  all  the  possible 
scattered  fields  can  be  included.  In  the  GTD  scattering  code  the 
importance  of  the  neglected  terms  are  determined  by  the  size  of  the 
so-called  gliches  or  jumps  ■'n  the  pattern  trace.  If  the  gliches 
are  small  no  additional  terms  are  needed  for  a  good  engineering  solu¬ 
tion.  If  the  gliches  are  large  if  may  be  necessary  to  include  more 
terms  in  the  solution.  In  any  case  the  user  has  a  gauge  with  which 
he  can  examine  the  accuracy  of  the  results  and  is  not  falsely  led 
into  believing  a  result  is  correct  when  in  fact  there  could  be  an 
error.  The  examples  in  Chapter  VI  illustrate  these  points  and  con¬ 
firm  the  validity  of  the  solution. 
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The  source  presently  considered  in  the  computer  code  is  an  elec¬ 
tric  or  maonetic  radiator  with  a  cosine  distribution  in  one  dimension 
and  either  a  zero  width  or  a  uniform  distribution  in  the  other.  It 
has  arbitrary  length  and  width,  arbitrary  magnitude  and  phase,  and 
arbitrary  orientation.  The  source  distribution  is  given  as  follows 


dipole  source: 


x  =0,  y  =0, 
P  P 


aperture  source: 


J(vV 

M(zp,xp) 


v°> 


-HAWS  HAWS 

2  -  xp  -  2  ’ 


-HS  ,  HS 

T^zpir 


where  the  x  dimension  is  oriented  in  the  THOX  and  PHOX  direction 
P 

and  the  zp  dimension  is  oriented  in  the  THOZ  and  PHOZ  direction, 
as  illustrated  in  Figure  5.  The  far-zone  electric  field  is  given 
by 


(t  )  =  E  F  (fl  )F  f P  ,4>  1 

p '  0  7  0  X  p 


where  for  an  electric  source 


9n  ^  I  HS. 

p  7T  m 


E  = 
o  i 


0  ^  J  HS  HAWS, 

^  p  n  m 


dipole  source 


aperture  source, 


’0 


and  for  a  magnetic  source 


'•'1 


-*p  J  «• 


-<t>  ^  M  HS  HAWS, 

p  7T  m 


dipole  source 


aperture  source, 


and  where 


sin6  cosfirHS  cose  ) 

fz«u  *  — <4 — ? - p- 

p  ( 1-4HS  COS0J 


%<VV  "1 


s  i  n  f  7T  HAWS  sinB  cos*  ) 

_  P  P 

HAWS  sinfl  cos* 

P  P 


,  dipole  source 


,  aperture  source. 


Any  arbitrary  antenna  can  he  simulated  by  superposition  of  the  ele¬ 
ments  hy  mak'r.q  th(  length  HS  small  f HS  =  OJa)  and  HAWS=0,  spacing 
the  elements  less  than  a  quarter  wavelength,  and  then  weighting  thei>" 
magnitudes  and  phases  to  simulate  the  current  distribution  of  the 
desired  antennaf9].  Since  the  radiation  pattern  is  relatively  in¬ 
sensitive  to  the  current  distribution  this  method  works  very  well. 
This  current  distribution  information  can  be  obtained  using  the  NEC- 
Moment  Method  Code.  Using  that  approach  the  field  from  an  individual 
element  in  the  moment  method  interface  section,  assuming  HS  is  small, 
is  given  by 


‘m 

2\ 


sin8 ' 


e-jkr* 

r5 


1 1 


s 


>>  •  >  < 1  s  consistent  resultant  fields  between  the  Moment  Method 

and  Pasic  Scattering  Cod^s.  It  should  be  emphasized  that  the  time 
required  to  calculate  a  radiation  pattern  increases  by  a  factor  N 
where  N  is  the  nurnbe*-  of  elemental  radiators  specified. 

Even  though  the  present  code  is  based  on  the  above  source  model 
it  can  be  easily  changed  by  modifying  the  SOURCE  subroutine  and  the 
SOURCP  subroutine,  which  contains  the  derivative  of  the  pattern  for 
the  slope  diffracted  fields.  This  information  is  given  in  more  de¬ 
tail  in  the  Code  Manual. 

The  brief  discussion  of  the  operation  of  the  scattering  code 
qiven  above  should  help  the  user  get  a  feel  for  the  overall  code 
so  he  might  better  understand  the  codes  capabilities  and  interpret 
its  results.  The  code  is  designed,  however,  so  that  the  general 
user  can  run  the  code  without  knowinq  all  the  details  of  its  operation. 
Yet,  he  must  become  familiar  with  the  input/output  details  which 
will  be  discussed  in  the  next  three  chapters. 
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CHAPTER  III 

DEFINITION  OF  INPUT  DATA 


The  method  used  to  input  data  into  the  computer  code  is  presently 
based  on  a  command  word  system.  This  is  especially  convenient  when 
more  than  one  problem  is  to  be  analyzed  during  a  computer  run.  The 
code  stores  the  previous  input  data  such  that  one  need  only  input 
that  data  which  needs  to  be  changed  from  the  previous  execution. 

Also,  there  is  a  default  list  of  data  so  for  any  given  problem  the 
amount  of  data  that  needs  to  be  input  has  been  shortened.  The  command 
word  options  presently  available  are  listed  in  Table  2. 


Command 


Table  2 

Input  Command  Options 


Description 


Page 


AM 

CE 

CG 

CM 

EN 

FR 

GP 

LP 

NC 

NG 

NP 

NS 

NX 

PD 

PG 

PP 

PR 

RG 

RT 

SG 

TO 

UN 

US 

XQ 


NEC  nr  AMP  Input  44 
Last  or  Only  Comment  Card  15 
Cylinder  Geometry  Input  36 
Comment  Card  15 
End  of  Execution  57 
Frequency  24 
Infinite  Ground  Plane  35 
Line  Printer  Listing  of  Results  54 
No  Cylinder  51 
No  Ground  Plane  50 
Next  Set  of  Plates  49 
Next  Set  of  Sources  52 
Next  Problem  53 
Pattern  Data  Desired  25 
Plate  Geonetry  Input  32 
Pen  Plot  Results  55 
Pov.*°r  Radiated  Input  48 
Far  Field  Range  Input  28 
Translate  and/or  Rotate  Coordinates  29 
Source  Geometry  Input  39 
Test  Data  Generation  Options  17 
Units  of  Input 

Units  of  HS  and  HAWS  in  SG:  ?? 
Execute  Program  56 


In  this  system,  all  linear  dimensions  may  be  specified  in  either 
in»'f>>rs,  inches,  feet  and  all  angular  dimensions  are  in  degrees. 

All  the  d  miens  ions  are  eventually  referred  to  a  fixed  cartesian  co¬ 
ordinate  system  used  as  a  common  reference  for  the  source  and  scat¬ 
tering  structures.  There  is,  however,  a  geometry  definition  coordi¬ 
nate  system  that  may  be  defined  using  the  “RT : "  command.  This  com¬ 
mand  enables  the  user  to  rotate  and  translate  the  coordinate  system 
to  be  used  to  input  any  selected  data  set  into  the  best  coordinate 
system  for  that  particular  geometry.  Once  the  "RT:"  command  is  used 
all  the  input  following  the  command  will  he  in  that  rotated  and  trans¬ 
lated  coordinate  system  until  the  "RT:"  command  is  called  again. 

See  below  for  more  details.  There  is  also  a  separate  coordinate 
system  that  can  be  used  to  define  a  pattern  coordinate  system.  This 
is  discussed  in  more  detail  helow  in  terms  of  the  "PD:"  command. 

It  is  felt  that  the  maximum  usefulness  of  the  computer  code 
can  be  achieved  using  it  on  an  interactive  computer  system.  As  a 
consequence,  all  input  data  are  defined  in  free  format  such  that 
the  operator  need  only  put  commas  between  the  various  inputs.  This 
allows  the  user  on  an  interactive  terminal  to  avoid  the  problems 
associated  with  typing  in  the  field  length  associated  with  a  fixed 
format.  This  method  also  is  useful  on  batch  processing  computers. 

Note  that  all  read  statements  are  made  on  unit  #5,  i.e.,  READ  (5,*), 
where  the  symbol  refers  to  free  format.  Other  machines,  however, 
may  have  different  symbols  representing  free  format. 

In  all  the  following  discussions  associated  with  logical  vari¬ 
ables  a  "T"  will  imply  true,  and  an  "F"  will  imply  false.  The  com¬ 
plete  words  true  and  false  need  not  be  input  since  most  compilers 
iust  consider  the  first  character  in  determining  the  state  of  the 
lonical  variable. 
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The  following  list  defines  in  detail  each  command  word  and  the 
variables  associated  with  them.  Chapter  VI  will  give  specific  examples 
using  this  input  method. 

A.  Commands  CM:  and  CE: 


These  commands  enable  the  user  to  place  comment  cards  in  the 
input  and  output  data  in  order  to  help  identify  the  computer  runs 
for  present  and  future  reference. 


1.  READ:  (IRfl) ,  1=1,24) 


a)  IR ( I ) :  This  is  an  integer  dimensioned  array  used  to 

store  the  command  word  and  comments.  Each  card 
should  have  CM:  or  CE:  on  them  followed  by  an 
alphanumeric  string  of  characters.  The  CM:  com¬ 
mand  implies  that  there  will  be  another  comment 
card  following  it.  The  last  comment  card  must 
have  the  CE:  command  on  it.  If  there  is  only 
one  comment  card  the  CE:  command  must  be  used. 

Note:  It  is  possible  to  place  comments  to  the  right  of  all  the  command 
words,  if  desired. 


FLOW  DIAGRAM  FOR  CE: 
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FLOW  DIAGRAM  FOR  CM: 


B.  Command  TO: 


This  command  enables  the  user  to  obtain  an  extended  output  of 
various  intermediate  quantities  in  the  computer  code.  This  is  use¬ 
ful  in  testing  the  program  or  in  analyzing  the  contributions  from 
various  scattering  mechanisms  in  terms  of  the  total  solution. 

1.  READ:  LDEBUG,  LTEST,  LOUT 

a)  LDEBUG:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  debug  the  program  if 
errors  are  suspected  within  the  program. 

If  set  true,  the  program  prints  out  data 
on  unit  #6  associated  with  each  of  its  in¬ 
ternal  operations.  These  data  can,  then 
he  compared  with  previous  data  which  are 
known  to  be  correct.  It  is,  also,  used 
to  insure  initial  operation  of  the  code. 

Only  one  pattern  angle  is  considered. 

( normally  set  false) 

b)  LTEST:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  test  the  input/output 
associated  with  each  subroutine.  The  data 
written  out  on  unit  #6  are  associated  with 
the  data  in  the  window  of  the  subroutine. 

They  are  written  out  each  time  the  subroutine 
is  called.  It  is,  also,  used  to  insure 
initial  operation  of  the  code.  Only  one 
pattern  angle  is  considered. 

(normally  set  false) 

c)  LOUT:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  output  data  on  unit 
#6  associated  with  the  main  program.  It 
too  is  used  to  initially  insure  proper  opera¬ 
tion.  It  can,  also,  be  used  to  examine 
the  various  components  of  the  pattern. 

This  is  expeciallv  useful  to  someone  in¬ 
terested  in  analyzing  which  scattering  center 
contributed  in  a  particular  direction. 

See  Table  3  for  a  list  of  the  fields  and 
their  identifiers. 

?.  READ:  L SLOPE,  LCORNR ,  LSOR 

a)  LSLOPE:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  tell  the  code  whether 
or  not  slope  diffraction  is  desired  during 
the  computation. 

(normally  set  true) 
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b)  L COR NR: 


This  is  -3  logical  variable  defined  by  T  or 
F.  It  is  used  to  tell  the  code  whether  or 
not  corner  diffraction  is  desired  during  the 
computation. 

(normal ly  set  true) 

c)  LSOR:  This  is  a  logical  variable  which  is  defined 

by  T  or  F.  It  is  used  to  specify  whether 
or  not  the  operator  wants  simply  the  antenna 
pattern  alone. 

(normal ly  set  false) 

3.  READ:  JMN(l),  JMX(l),  JMN{2),  JMX(2),  JMN(3),  JMX(3) 

a)  JMN(l),  JMX(l):  These  are  integer  variables  used  to 

specify  a  set  of  individual  scattering  compon¬ 
ents  that  are  to  be  included  in  the  scattered 
field  computation  for  the  plate  structures 
alone.  JMN(l)  is  the  minimum  component  number 
and  JMX ( 1 )  is  the  maximum  component  number 
for  the  range  of  the  set  where  the  components 
are  defined  by  the  following  number  designations: 
0  =  skip  the  plates  section 

1  =  incident  field 

2  =  single  reflected  fields 

3  =  double  reflected  fields 

4  =  single  diffracted  fields 

5  =  ref lected-diffracted  fields 

6  =  diffracted-ref lected  fields 

7  =  identifies  double  diffracted  problem  areas 

(double  diffracted  fields  are  not  computed 
at  present). 

Normally  0MN(1)=1  and  JMX(1)=7.  This  would 
compute  all  the  available  field  values  for 
a  convex  or  concave  plate  structure. 

b)  JMN(2),  JMX(2):  These  are  integer  variables  used  to 

specify  a  set  of  individual  scattering  com¬ 
ponents  that  are  to  be  incluaed  in  the  scat¬ 
tered  field  computation  for  the  finite  elliptic 
cylinder  structure  alone.  JMN(2)  is  the 
minimum  component  number  and  JMX ( 2 )  is  the 
maximum  component  number  for  the  range  of 
the  set  where  the  components  are  defined  by 
the  following  number  designations: 

0  =  skip  the  cylinder  section 

1  =  incident,  reflected,  transition  and  creeping 

wave  fields. 

2  =  single  reflected  fields  from  the  end  caps. 
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3  =  single  diffracted  fields  from  the  end 
cap  rims. 

Normally  JMN(2)=1  and  JMX(2)=3.  This  would 
compute  all  the  field  values  for  a  finite 
elliptic  cylinder  structure. 

c)  JMN(3),  JMX(3):  These  are  integer  variables  used  to 

specify  a  set  of  individual  scattering  compon¬ 
ents  that  are  to  be  included  in  the  scattered 
field  computation  for  the  interactions  between 
the  plate  and  cylinder  structures.  JMN(3) 
is  the  minimum  component  number  and  JMX (3) 
is  the  maximum  component  number  for  the  range 
of  the  set  where  the  components  are  defined 
by  the  following  number  designations: 

0  =  skip  the  plate-cylinder  interaction  section. 

1  =  fields  reflected  from  the  plates  then 

reflected  or  diffracted  from  the  cylinder. 

2  =  fields  reflected  or  diffracted  from  the 

cylinder  then  reflected  from  the  plates. 

3  =  fields  reflected  from  the  cylinder  then 

diffracted  from  the  plates. 

4  =  fields  diffracted  from  the  plates  then 

reflected  from  the  cylinder. 

Normally  JMN(3)=1  and  JMX(3)=4. 

FLOW  DIAGRAM  FOR  TO: 


READ:  LDEBUG,LTEST,LOUT 


WRITE:  L DEBUG, LTEST, LOUT 


READ:  LSLOPE.LCORNR.LSOR 


WRITE:  LSLOPE ,LC0RNR ,LSOR 


READ:  JMN ( 1 ) , JMX ( 1 ) , JMN ( 2 ) , JMX ( 2 ) , JMN ( 3 ) , JMX ( 3 ) 


WRITE:  JMN { 1 ) , JMX ( I ) , JMN ( 2 ) , JMX ( 2 ) , JMN ( 3 ) , JMX ( 3 ) 


Table  3 

Individual  Field  Types  Printed  when  LOUT=.TRUE 


L 

K 

J 

r 

Field  Type 

100 

0 

0 

0 

Direct  field  when  plates  are  present 

200 

MP 

0 

0 

Field  reflected  from  plate  MP 

300 

MP 

MPP 

0 

Field  reflected  from  plate  MP  then 
reflected  from  plate  MPP 

600 

MP 

ME 

0 

Field  diffracted  from  edge  ME  of 
plate  MP 

650 

MP 

ME 

0 

Field  diffracted  frcm  the  corners 
of  edge  ME  of  plat?  mp 

700 

MR 

MP 

ME 

Field  reflected  from  plate  MR  then 
diffracted  from  edge  ME  of  plate  MP 

750 

MR 

MP 

ME 

Field  reflected  frcm  plate  MR  then 
diffracted  by  the  corners  of  edge 

ME  of  plate  MP 

300 

MP 

ME 

MR 

Field  diffracted  from  edge  ME  of 
plate  MP  then  reflected  f-om  plate 

Field  diffracted  from  the  corners 
of  edge  me  of  plate  mp  then  reflected 
from  plate  MR 

350 

MP 

ME 

MR 

110 

3 

0 

0 

Direct  field  ?hen  only  cylinders 
alone  are  present 

120 

0 

0 

0 

Geometrical  optics  field  reflected 
by  cylinder  (for  comparison  only) 

130 

0 

0 

0 

Field  scattered  by  the  curved  sur¬ 
face  of  the  cylinder 

150 

MC 

0 

0 

Field  reflected  by  end  cap  MC  of 
the  cylinder 

500 

MC 

0 

0 

Field  diffracted  by  the  end  cap 
rim  MC  of  the  cylinder 

240 

VP 

0 

0 

Geometrical  optics  f:eld  reflected 
from  plate  MP  then  reflected  from 
the  curved  surface  of  the  cylinder. 
(For  comparison  only) 

250 

MP 

0 

0 

Field  reflected  from  plate  MP  and 
then  scattered  by  the  curved  sur¬ 
face  of  the  cyl inder 

410 

MP 

0 

0 

Geometrical  optics  field  reflected 
from  the  curved  surface  of  the  cy¬ 
linder  and  then  reflected  from  plate 
MP,  (For  comparison  only) 

420 

MP 

0 

0 

Field  scattered  from  the  curved 
surface  of  the  cylinder  then  re¬ 
flected  from  olate  MP 

?40 

MP 

ME 

0 

Field  reflected  f-m  the  curved 
surface  of  the  cylinder  then  dif¬ 
fracted  by  edge  ME  of  plate  MP 

950 

MP 

ME 

0 

Field  diffracted  from  edge  ME  of 
p’ate  MP  then  reflected  from  the 
curved  surface  of  the  cylinder 

I A\GL E 

1  ANGLE 

INDEX 

INDEX 

Sum  of  fields  of  a  given  type  (INDEX) 
for  a  given  angle  ' [ANGLE) 

1000 

1 ANGLE 

I  ANGLE 

1  ANGLE 

Total  field  for  a  given  angle  (JANGLE) 

?0 
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C.  Command  UN: 


This  command  enables  the  user  to  specify  the  units  of  all  the 
linear  dimensions  to  be  input  after  the  command  is  called.  (The  one 
exception  is  the  source  length  HS  and  width  HAWS,  see  command  US:) 


READ:  I UN I T 
a)  IUNIT: 

IUNIT  = 


This  is  an  integer  variable  that  indicates 
the  units  for  the  input  data  that  follows, 
such  that  if 
[  1  -►  meters 
l  2  ■*  feet 
1  3  -►  inches 


FLOW  DIAGRAM  FOR  UN: 


READ:  IUNIT 


Set  up  Conversion  constant 

Write  units  specified 


0.  Command  US: 


This  command  enables  the  user  to  specify  the  units  of  the  source 
>ngth  HS  and  width  HAWS  to  be  input  after  the  command  is  called.  These 
variables  are  in  the  command  SG: . 


1.  READ:  IUNST 


a)  IUNST: 


IUNST  • 


This  is  an  integer  variable  that  indicates 
the  units  for  the  input  data  HS  and  HAWS  that 
follows,  such  that  if 
[0  ■*  wavelengths 
/  1  ■*  meters 
|  2  -*■  feet 
3  -►  inches 


Note  that  if  the  units  are  specified  to  be  wavelengths  for  one  source 
it  must  be  wavelengths  for  all  of  the  sources  specified. 


! 

I 
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FLOW  DIAGRAM  FOR  US: 


I 

I 

I 
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Command  FR: 


This  command  enables  the  user  to  define  the  frequency  in  gigahert 
1.  READ:  FRQG 

a)  FRQG:  This  is  a  real  variable  which  is  used  to 

define  the  frequency  in  gigahertz. 

FLOW  DIAGRAM  FOR  FR: 


F. 


Command  PD: 


This  command  enables  the  user  to  define  the  pattern  coordinate 
system,  the  pattern  cut,  and  the  angular  range  that  is  desired.  The 
geometry  is  illustrated  in  Figure  1. 


1.  READ:  THCZ,  PHCZ,  THCX,  PHCX 

a)  THCZ, PHCZ:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  z  -axis  of  the  pattern  coordinate  system 
as  ifpit  was  a  rad'al  vector  in  the  reference 
coordinate  system. 

b)  THCX, PHCX:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  x-axis  of  the  pattern  coordinate  system 
as  ifpit  was  a  radial  vector  in  the  reference 
coordinate  system. 


Note  that  the  new  xp-axis  and  zp-axis  must  be  defined  orthogonal  to 
each  other.  The  new  yp-axis  is  found  from  the  cross  product  of  the 
xp-  and  zp-axis. 


2.  READ:  LCNPAT,  TPPO 


a)  LCNPAT: 
LCNPAT  = 


This  is  a  logical  variable  that  defines  the 
pattern  cut  desired,  such  that  if 
>T  -  THETA  CUT (conic  cut) 

\F  -  PHI  CUT  (PHI  constant) 


b)  TPPD: 

LCNPAT  = 


This  is  a  real  variable  that  defines  the 
pattern  angle  that  is  to  be  held  constant, 
such  that  if 

/ T  ■*  TPPD  =  THP  constant 
\F  •*  TPPD  =  PHP  constant 


3.  READ:  IB,  IE,  IS 


a)  18, IE, IS:  These  are  integer  variables  used  to  define 
angles  in  degrees.  They  are,  respectively, 
the  beginning,  ending,  and  incremental  values 
of  the  pattern  angle. 


FLOW  DIAGRAM  FOR  PD: 


THCZ,PHCZ,THCX,PHCX 


Define  the  pattern  cut  x  and  z  axes 
from  the  input  data  coordinate  system 


^  Are  the.yNxSi 
perpendicular 
to  each  othej^ 


Write  warning  message 


Define  the  y  axis  from  a  cross 
product  of  the  x  and  z  axes 


Write  the  x,y,z  axes  of  the  pattern 
cut  coordinate  system 


READ:  ICNPAT.TPPD 


Write  the  pattern  cut  desired 


READ:  IB, IE, IS 


WRITE:  IB, IE, IS 
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i 
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Figure  !a.  Definition  of  pattern  coordinate  system. 


XP 


Figure  lb.  Conic  pattern  cut,  LCNPAT=.TRUE. ,  tppd=THP. 


*P 

Figure  1c.  Constant  Phi  pattern  cut,  LCNPAT=. FALSE. ,  TPPD=PHP. 
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Command  RG : 


This  command  enables  the  user  to  specify  a  far  field  distance, 

R,  to  the  observer.  The  fields  are  then  normalized  by  the  factor 
exp(-jkR)/R. 

1.  READ:  RANGS 

a)  This  is  a  real  variable  which  is  used  to  specify  the 
far  field  range,  R. 

Note  that  R  should  be  in  the  far  field  of  the  scattering  structure, 
that  is,  R  :■  2D  />.  where  D  is  the  maximum  dimension  of  the  structure. 
However,  if  R  _  10^,  then  the  factor  exp(-jkR)/R  is  suppressed. 

FLOW  DIAGRAM  FOR  RG: 


T 


H.  Command  RT: 


This  command  enables  the  user  to  translate  and/or  rotate  the  coordin 
ate  system  used  to  define  the  input  data  in  order  to  simplify  the  speci¬ 
fication  of  the  plate,  cylinder,  and  source  geometries.  The  geometry 
is  illustrated  in  Figure  2. 


1.  READ:  ( TR ( N ) ,  N=1 ,3) 

a)  TR ( N ) :  This  is  a  dimensional  real  variable.  It  is 

used  to  specify  the  origin  of  the  new  coordin¬ 
ate  system  to  be  used  to  input  the  data  for 
the  source  or  the  scattering  structures. 

It  is  input  on  a  single  line  with  the  real 
numbers  being  the  x,y,z  coordinates  of  the 
new  origin  which  corresponds  to  N-1,2,3, 
respectively. 

2.  READ:  THZP,  PHZP,  THXP,  PHXP 

a)  THZP,  PHZP:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  z-axis  of  the  new  coordinate  system  as 
if  it  was  a  radial  vector  in  the  reference 
coordinate  system. 

b)  THXP,  PHXP:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  x-axis  of  the  new  coordinate  system  as 
if  it  was  a  radial  vector  in  the  reference 
coordinate  system. 

The  new  x-axis  and  z-axis  must  be  defined  orthogonal  to  each  other. 
The  new  y-axis  is  found  from  the  cross  product  of  the  x-  and  z-axes. 

All  the  subsequent  inputs  will  be  made  relative  to  this  new  coordinate 
system,  which  is  shown  as  x^,  yt,  zt  ,  unless  command  "RT:"  is  called 
again  and  redefined. 
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FLOW  DIAGRAM  FOR  RT: 


Define  the  y  axis  from  across 
product  of  the  x  and  z  axes 


30 


I. 


Command  PG: 


This  command  enables  the  user  to  ut. ine  the  geometry  of  the  flat 
plate  structures  to  be  considered.  The  geometry  is  illustrated  in 
Figure  3.  One  call  to  this  command  defines  one  plate.  The  number  of 
plates  in  the  structure  are  automatically  counted  by  the  number  of  calls 
to  this  command. 


1.  RE AO:  MEP(MP) 

a)  MEP  (MP):  This  is  a  dimensioned  integer  variable.  It 
is  used  to  define,  the  number  of  corners  (or 
edges)  on  the  MPtn  plate. 

2.  REAO:  (XX(MP,ME,N) ,  N=l,3) 

a)  XX(MP,M£,N):  This  is  a  triply  dimensioned  real  variable. 

It .Is  used  to  specifyfthe  location  of  the 
MEtn  corner  of  the  MPtn  plate.  It  is  input 
on  a  single  line  with  the  real  numbers  being 
the  x,y,z  coordinates  of  the  corner,  in  the 
specified  coordinate  system,  which  correspond 
to  N=1 ,2,3,  respectively,  in  the  array.  For 
example,  the  array  will  contain  the  following 
for  plate  #1  and  corner  #2  located  at  x-2., 
y=4. ,  z=6. : 

xxn,2,n=p. 

xxn,?,?M. 

XXM,2,3)«fi. 

This  data  is  input  as:  2., 4. ,6. 


This  read  statement  will  be  called  MEP(MP)  times  so  that  all  the 
corners  are  defined.  As  an  example,  the  input  data  for  the  flat  plate 
structure  given  in  Figure  3,  is  given  by 


4 

1.,  1.,  0 
-1 . ,  1 . ,  0. 
-1 . ,-l . ,  0. 
1 . ,-l . ,  0. 


:number  of  corners  for  plate  #1 
:  corner  #1 
:  corner  #2 
:  corner  #3 
:  corner  #4  • 


See  Chapter  IV  for  further  details  on  how  to  number  the  corners. 


32 


Note  that  the  program  will  keep  increasing  the  number  of  plates  in  the 
solution  by  the  number  of  calls  to  this  command  unless  the  NP:  or  NX: 
commands  are  called  to  reinitialize  the  plate  geometry. 


Presently:  1  1  MP  1  14 

1  £  ME  £  6 

1  <  N  <  3 


Figure  3.  Definition  of  flat  plate  geometry. 
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>1.  Command  GP : 


This  command  enables  the  user  to  specify  a  perfectly-conducting 
infinite  ground  plane  in  the  x^-y^  plane. 


FLOW  DIAGRAM  FOR  GP: 


***  . 


K.  Command  CG: 

This  command  enables  the  user  to  define  the  geometry  of  the  finite 
elliptic  cylinder  structure  to  be  considered.  Note  only  one  may  be 
specified.  The  geometry  is  illustrated  in  Figure  4. 


READ:  AA,  BB 
a)  AA: 


b)  BB: 


This  is  a  real  variable  which  defines  the 
radius  of  the  elliptic  cylinder  on  the  x^-axis 
of  the  cylinder. 

This  is  a  real  variable  which  defines  the 
radius  of  the  elliptic  cylinder  on  the  yt~axis 
of  the  cylinder. 


READ:  ZCN,  THTN,  ZCP,  THTP 

a)  ZCN:  This  is  a  real  variable  that  defines  the 

position  of  the  center  of  the  most  negative 
endcap  on  the  z^-axis  of  the  cylinder. 

b)  THTN:  This  is  a  real  variable.  It  is  input  in  de¬ 

grees  and  defines  the  angle  the  surface  of 
the  most  negative  endcap  makes  with  the  nega¬ 
tive  zt-axis  in  the  xt-zt  plane. 

c)  ZCP:  This  is  a  real  variable  that  defines  the 

position  of  the  center  of  the  most  positive 
endcap  on  the  z^-axis  of  the  cylinder. 

d)  THTP:  This  is  a  real  variable.  It  is  input  in  degrees 

and  defines  the  angle  the  surface  of  the  most 
positive  endcap  makes  with  the  positive  zt-axis 
in  the  x  -zt  plane. 
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j 


r 


Figure  4.  Definition  of  finite  elliptic  cylinder  geometry. 
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Command  SG: 


This  command  enables  the  user  to  spec’fv  the  location  and  type 
of  source  to  be  used.  The  geometry  is  illustrated  in  Figure  5.  One 
call  to  this  command  defines  one  source.  The  number  of  sources  in  the 
problem  are  automatically  counted  by  the  number  of  calls  to  this  command. 

1.  READ:  (XSS  (MS,N)  N=l,3) 

a)  XSS  { MS , N ) :  This  is  a  doubly  dimensioned  real  array 

which  is  used  to  define  the  x,y,z  location 
of  the  MScn  elemental  radiator  in  the  speci¬ 
fied  cartesian  coordinate  system.  Again, 
a  single  line  of  data  contains  the  x,y,z 
( N= 1,2,3)  locations. 

2.  REAO:  THOZ(MS),  PHOZ(MS),  THOX(MS),  PHOX(MS) 

a)  TH0Z(MS),PH0Z(MS) :  These  are  real  arrays  which  are 

used  to  define  the  orientation  of  the  MStn 
element  in  the  specified  cartesian  coordinate 
system.  They  are  input  in  degrees,  as  spherical 
angles,  that  define  a  radial  direction  which 
is  parallel  to  the  MStn  element  current  flow 
for  a  dipole  antenna  or  which  is  parallel 
to  the  length  of  an  aperture  antenna. 

b)  TH0X(MS),PH0X(MS):  These  are  real  arrays  which  are. 

used  to  define  the  orientation  of  the  MScn 
element  in  the  specified  cartesian  coordinate 
system.  They  are  input  in  degrees,  as  spherical 
angles,  that  define  a.radial  direction  which 
is  parallel  to  the  MStn  elements  aperture 
width  or  which  is  parallel  to  a  slots  width. 

For  a  dipole  antenna,  these  angles  can  be 
made  in  a  convenient  direction. 

The  x-axis  and  z-axis  specified  by  these  angles  must  be  defined 
orthogonal  to  each  other.  The  y-axis  is  found  by  the  cross  product 
of  the  x-  and  z-axes. 
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3.  READ:  IMS(MS),  HS(MS),  HAWS(MS) 


a)  IMS(MS) :  This  is  an  integer  array  which  is  used  to 

define  whether  the  MS™  source  is  an  electric 
or  magnetic  elemental  radiator. 

IMS(MS)  =  0  -  electric 
IMS(MS)  =  1  -*  magnetic 

b)  HS(MS):  This  is  a  real  array. which  is  used  to  input 

the  length  of  the  MS™  element. 


c)  HAWS (MS)  This  is  a  real  array.which  is  used  to  input 
the  width  of  the  MS™  element  in  the  case 
of  an  aperture  antenna.  If  HAWS(MS)=0  then 
it  is  assumed  to  be  a  dipole  antenna. 


Note  that  the  units  of  the  variables  HS(MS)  and  HAWS(MS)  can  be 
specified  by  the  US:  command.  If  wavelength  is  chosen  as  the  units 
then  all  the  sources  must  be  specified  in  wavelengths. 


4.  READ:  WM(MS),  WP(MS) 

a)  WM(MS),  WP(MS):  These  are  real  dimensioned  arrays  used 
to.define  the  excitation  associated  with  the 
MS™  element.  The  magnitude  is  given  by  WM 
and  the  phase  in  degrees  by  WP. 


Note  that  the  program  will  keep  increasing  the  number  of  sources 
in  the  solution  by  the  number  of  calls  to  this  command  unless  the  NS: 
or  NX:  commands  are  called  to  reinitialize  the  source  geometry. 

Presently,  1  <  MS  £  50 
1  <  N  <  3 
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^  If  ' 
IUNST=0 


WRITE:  HS(MS) ,HAWS(MS) 


Convert  source  height  and  length 
to  meters 


WRITE:  HS'MS).HAWS(MS) 


/  WRITE:  WM(MS).WPfMS) 


Define  the  y  axis  from  across 
product  of  the  x  and  z  axes 


WRITE  x,y,z  axes  of 
source  coordinate  system 


PHOZ 


M.  Command  AM: 


This  command  enables  the  user  to  interface  the  Numerical  Electro¬ 
magnetics  Code  (NEC)-Moment  Method  Code  with  the  Basic  Scattering 
Code  in  order  to  use  the  antenna  modeling  capabilities  of  NEC  to 
specify  the  needed  source  data  such  as  location  and  current  weights 
of  the  elements.  The  geometry  is  illustrated  in  Figure  f>. 


1.  READ:  PRAD 

a)  PRAD:  This  is  a  real  variable  which  is  used  to  define 

the  total  power  radiated  in  watts  from  the 
NEC  modeled  antenna.  This  allows  the  directive 
gain  to  be  computed.  If  P.  is  substituted 
for  Pra(j  the  power  gain  wilT  be  computed. 

2.  READ:  MSX 

a)  MSX:  This  is  an  integer  variable  which  defines 

the  maximum  number  of  elemental  wire  radiators 
that  have  been  used  in  the  NEC  code  to  model 
the  antenna. 

3.  READ:  (XSS(M$,N),N=1,3),  HS(MS),  THOZ(MS),  PHOZ(MS) 

a)  XSS(MS.N):  This  is  a  doubly  dimensioned  real  array  which 

istused  to  define  the  x,y,z  location  of  the 
MStn  elemental  radiator  in  the  specified 
cartesian  coordinate  system. 

b)  HS(MS):  This  is  a  real  array. which  is  used  to  input 

the  length  of  the  MStn  element  in  the  units 
specified  by  I UN IT  from  the  UN:  command  or 
from  the  default  option. 

c)  THOZ(MS),  PHOZ(MS):  These  are  real  arrays  which  are 

used  to  define  the  orientation  of  the  MStn 
element  in  the  specified  cartesian  coordinate 
system.  The  THOZ,  PHOZ  angles  are  in  degrees 
and  define  a  radi aT  direction  which  is  paral¬ 
lel  to  the  MSin  element  current  flow.  The 
angle  THO  is  the  angle  of  the  element  measured 
up  from  the  x-y  plane.  The  angle  PHO  is  the 
phi  angle  in  a  normal  spherical  coordinate 
system. 


4.  READ:  WM(MS),  WP(MS) 

a)  WH(MS),  WP(MS):  There  are  real  dimensioned  arrays  used 
totdefine  the  excitation  associated  with  the 
MSin  element.  The  real  part  is  given  by  WM 
and  the  imaginary  part  by  WP. 

Presently,  1  <  MS  ±  50. 

1  <  N  1  3 


Note  that  for  the  NEC  code  input  all  the  elements  are  assumed 
to  be  electric  current  elements. 


Figure  5.  Definition  of  NEC  geometry. 


FLOW  DIAGRAM  FOR  AM: 


READ:  (WM(MS) ,WP(MS) ,MS=1 ,MSX) 


|  Rotate  and  translate  coordinates 


Change  units  to  meters 


Define  NEC  source 
coordinate  systems 


L 


Write  new  and  old 
coordinates  and  units 
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N. 


Command  PR: 


This  command  enables  the  user  to  specify  the  total  power  radiated 
by  the  antenna  or  the  input  power  to  the  antenna. 

1.  READ:  PRAD 

a)  PRAO:  This  is  a  real  variable.  It  is  input  in  watts 

and  defines  the  total  power  radiated  by  the 
antenna  (or  input  power  to  the  antenna). 

-in 

Note  that  if  PRAD  _<  10  ,  it  will  be  assumed  that  the  power  radi¬ 

ated  (or  input  power)  was  not  specified. 


FLOW  DIAGRAM  FOR  PR: 


0.  Command  NP: 


This  command  enables  the  user  to  initialize  the  plate  data. 
All  of  the  plates  are  removed  from  the  problem  unless  they  are  re 
specified  following  this  command. 

FLOW  DIAGRAM  FOR  NP: 
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P.  Command  NG: 


This  command  enables  the  user  to  initialize  the  infinite  ground 
plane.  The  ground  plane  is  removed  from  the  problem  unless  it  is 
respecified  following  this  command. 


FLOW  DIAGRAM  FOR  NG: 


Set  ground  plane  logic  false 
LGRND=. FALSE. 


/  Write  message  / 


Q.  Command  NC : 


This  command  enables  the  user  to  initialize  the  cylinder  data. 
The  cylinder  is  removed  from  the  problem  unless  it  is  respecified 
following  this  command. 

FLOW  DIAGRAM  FOR  NC: 
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R.  Command  NS: 


This  command  enables  the  user  to  initialize  the  source  data. 
All  of  the  sources  are  removed  from  the  problem  unless  they  are  re 
specified  following  this  command. 


FLOW  DIAGRAM  FOR  NS: 


© 


_ 

Set  source 
LAMP=. 

_ 

logic  false 
FALSE. 

Reset  sour 
MS 

ce  counter 

X=0 

■ 

/ - - - 7 


Write  message 


S.  Command  NX: 


r« 

1 

1 

I  This  command  enables  the  user  to  reinitialize  the  commands  to 

■  |  their  default  conditions  specified  in  the  list  at  the  beginning  of 

[  the  main  program. 


T.  Command  LP: 


This  command  enables  the  user  to  specify  whether  a  line  printer 
listing  of  the  results  is  desired. 


1.  READ:  L WRITE 

a)  LWRITE:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  indicate  if  a  line 
printer  listing  of  the  total  fields  (Eft  , 

E,  )  is  desired.  p 

<i>P 

(normally  set  true) 

FLOW  DIAGRAM  FOR  LP: 
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U.  Command  PP: 


This  command  enables  the  user  to  specify  whether  a  polar  plot 
of  the  results  are  desired. 

1.  READ:  LPLT 


a)  LPLT:  This  is  a  logical  variable  defined  by  T 

or  F.  It  is  used  to  indicate  if  a  polar 
plot  of  the  total  fields  (Efl  ,EA  )  is  de¬ 
sired.  If  LPLT  is  false  thepres°  of  the 
READ  statements  for  this  command  will  be 
skipped. 

2.  READ:  RADIUS, IPLT 


a)  RADIUS:  This  is  a  real  variable  that  is  used  to 

specify  the  radius  of  the  polar  plot. 

b)  IPLT:  This  is  an  integer  variable  that  indicates 

the  type  of  polar  plot  desired,  such  that 
if 

Tp,  r-  fl  field  plot 
\2  ■*  power  plot 


j  c.  |;un i 

[3  -*•  dB  plot. 


The  fields  will  be  normalized  by  their  maximum  field  values. 


FLOW  DIAGRAM  FOR  PP: 
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This  command  is  used  to  execute  the  scattering  code  so  that 
the  total  fields  may  be  computed.  After  execution  the  code  returns 
for  another  possible  command  word. 


FLOW  DIAGRAM  FOR  XQ: 


/ 


n 


Write  message 


W.  Command  EN: 


This  command  enables  the  user  to  terminate  the  execution  of 
the  scattering  code. 


FLOW  DIAGRAM  FOR  EN: 


This  concludes  the  definition  of  all  the  input  parameters  to 
the  program.  The  program  would,  then,  run  the  desired  data  and  output 
the  results  on  unit  #6.  However  as  with  any  sophisticated  program, 
the  definition  of  the  input  data  is  not  sufficient  for  one  to  fully 
understand  the  operation  of  the  code.  In  order  to  overcome  this 
difficulty  the  next  chapter  discusses  how  the  input  data  are  inter¬ 
preted  and  used  in  the  program. 
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CHAPTER  IV 

INTERPRETATION  OF  INPUT  DATA 


This  computer  code  is  written  to  require  a  mimumum  amount  of 
user  information  such  that  the  burden  associated  with  a  complex  geometry 
will  be  organized  internal  to  the  computer  code,  For  example,  the 
operator  need  not  instruct  the  code  that  two  plates  are  attached 
to  form  a  convex  or  concave  structure.  The  code  flags  this  situation 
hy  recognizing  that  two  plates  have  a  common  set  of  corners  (i.e., 
a  common  edge).  So  if  the  operator  wishes  to  attach  two  plates  to¬ 
gether  he  needs  only  define  the  two  plates  as  though  they  were  isolated. 
However,  the  two  plates  will  have  two  identical  corners.  All  the 
geometry  information  associated  with  plates  with  common  edges  is 
then  generated  by  the  code.  The  present  code  also  will  allow  a  plate 
to  intersect  another  plate  as  shown  in  Figure  7.  It  is  necessary 
that  the  corners  defining  the  attachment  be  positioned  a  small  amount 
(approximately  10'^  wavelengths)  through  the  plate  surface  to  which 
it  is  being  connected. 

In  defining  the  plate  corners  it  is  necessary  to  be  aware  of 
a  subtlety  associated  with  simulating  convex  or  concave  structures 
in  which  two  or  more  plates  are  used  in  the  computation.  This  prob¬ 
lem  results  in  that  each  plate  has  two  sides.  If  the  plates  are 
used  to  simulate  a  closed  or  semi-closed  structure,  then  possibly 
only  one  side  of  the  plate  will  be  illuminated  by  the  antenna.  Con¬ 
sequently,  the  operator  must  define  the  data  in  such  a  way  that  the 
code  can  infer  which  side  of  the  plate  is  illuminated  by  the  antenna. 
This  is  accomplished  by  defining  the  plate  according  to  the  right- 
hand  rule.  As  one's  fingers  of  the  right  hand  follow  the  edges  of 
the  plate  around  in  the  order  of  their  definition,  his  thumb  should 
point  toward  the  illuminated  region  above  the  plate.  To  illustrate 
this  constraint  associated  with  data  format,  let  us  consider  the 
definition  of  a  rectangular  box.  In  this  case,  all  the  plates  of 
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Figure  7.  Data  format  used  to  define  a  flat  plate  intersecting 
another  flat  plate. 


the  box  must  be  specified  such  that  they  satisfy  the  right-hand  rule 
with  the  thumb  pointing  outward  as  illustrated  in  Figure  8.  If  this 
rule  were  not  satisfied  for  a  given  plate,  then  the  code  would  assume 
that  the  antenna  is  within  the  box  as  far  as  the  scattering  from 
that  plate  is  concerned. 

Another  situation  which  must  be  kept  in  mind  is  associated  with 
antenna  elements  mounted  on  a  plate.  The  code  automatically  deter¬ 
mines  that  the  antenna  element  is  mounted  on  the  plate.  It  assumes 
that  the  element  will  radiate  on  the  side  of  the  plate  into  which 

the  normal  points.  This  is  accomplished  in  the  code  by  automatically 

-5 

positioning  the  source  a  small  distance  (10  wavelengths  above  the 
plate  in  the  direction  of  the  normal  as  illustrated  in  Figure  9. 

It  is  important,  therefore,  to  follow  the  simple  rules  above  for 
defining  the  plate  normals  when  dealing  with  plate  mounted  antennas. 
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A 

n 


MONOPOLE 


EQUIVALENT 
DIPOLE  , 
LENGTH 

Li  c  .. 

8  <  I0_6X 


SLOT 

(MAGNETIC 
DIPOLE  ) 


“I- 

8  <  io_6X 


(o)  MONOPOLE  (b)  SLOT 


Figure  9.  Illustration  of  geometry  for  plate-mounted  antennas. 
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There  is,  also,  another  point  associated  with  antennas  mounted  on 
perfectly-conducting  flat  structures.  If  a  plate-mounted  monopole 
is  considered  in  the  computation,  one  should  input  the  equivalent 
dipole  length  and  not  the  monopole  length  (i.e.,  the  monopole  plus 
image  length  should  be  used  as  shown  in  Figure  9a).  The  code  auto¬ 
matically  handles  the  half  dipole  modes  associated  with  the  mono¬ 
pole.  The  plate-mounted  slot  is,  also,  automatically  taken  care 
of  by  the  code  as  shown  in  Figure  9b  if  a  magnetic  dipole  is  used. 

The  same  situation  arises  when  the  antenna  is  mounted  on  the 
elliptic  cylinder's  end  caps.  It  should  also  be  remembered  that 
the  antenna  can  not  be  mounted  on  the  curved  part  of  the  cylinder. 

In  general,  the  antenna  should  be  kept  a  wavelength  away,  however, 
this  can  often  be  relaxed  to  approximately  a  quarter-wavelength. 

In  the  present  code,  the  attachment  of  a  plate  corner  to  the 
curved  surface  of  the  cylinder  is  automatically  detected,  however, 
a  diffracted  field  from  the  plate-cylinder  junction  is  not  considered 
in  this  version.  If  the  plate-cylinder  junction  forms  a  straight, 
orthogonal  edge,  as  shown  in  the  aircraft  models  of  Figure  27,  image 
theory  alone  will  give  the  correct  results.  The  diffracted  fields, 
therefore,  are  not  needed.  If  the  plate-cylinder  junction  forms 
a  curved  edge  or  one  in  which  the  plate  and  cylinder  surface  are 
not  orthogonal  a  diffracted  field  from  that  edge  will  be  required 
in  the  solution.  This  will  be  added  when  time  and  effort  permit. 

Chapter  VI  has  a  set  of  sample  problems  to  illustrate  how  the 
operator  can  realize  the  versatility  of  the  code  and  still  satisfy 
the  few  constraints  associated  with  the  input  data  format. 


CHAPTER  V 

INTERPRETATION  OF  OUTPUT 

The  basic  output  from  the  computer  code  is  a  line  printer  listing 
of  the  results.  The  results  are  referenced  to  the  pattern  coordi¬ 
nate  system  that  was  described  in  Chapter  III  and  is  illustrated 
in  Figure  1.  Thus  the  total  electric  field  is  given  by 


*<VV 


Vep  *  *pE*p' 


The  fields  are  assumed  to  be  peak  values  given  in  volts/unit  when 

the  factor  e’^kR/R  is  suppressed  in  the  far  field.  If  an  R  value 

is  specified  using  the  "RG:"  command  then  the  results  will  be  in 

volts/meter.  The  results  are  displayed  in  three  sections.  The  first 

output  associated  with  the  E0p  field,  the  second  is  output  associated 

with  the  E  field,  and  the  third  is  output  associated  with  the  total 
<pP 

field.  The  first  section  is  displayed  as  follows:  0p,  4>p,  E0p, 

/EQp  (phase  of  E0p),  |EQp|  (magnitude  of  E0p),  Gd0  (directive  gain 

V-  lE9pl/|EepU»x  (nor,na'1zed  Eep),  Gd  nonn  (nor- 

malized  value  of  the  directive  gain).  The  second  section  is  similarly 
done  for  the  E^  field.  The  third  section  is  displayed  as  follows: 

B.  <t>  ,  G.  (directive  gain  of  E„4.J,  G  .  ^  (directive 

gain  of  Eminor)»  Y  (tilt  angle  of  polarization  ellipse,  axial  ratio. 


Gd  (total  directive  gain),  Grf  norm 
The  above  quantities  are  defined  as  follows: 


(normalized  total  directive  gain). 
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d  major 
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axial  ratio  = 


^minor 

^major 


lE^pl^ln^sin^I1 
I  E  pl  2  s1n?*coszY] 


n  =  free  space  impedance. 


The  value  P  ^  is  the  power  radiated.  It  is  input  into  the 
basic  scattering  code  through  the  AM:  command  for  the  NEC  moment 
method  data  or  through  the  PR:  command.  If  a  value  for  Pfacj  is  not 
given  to  the  code,  the  output  will  be  given  in  terms  of  the  radiation 
intensity  rather  than  the  directive  gain.  The  radiation  intensity, 

U,  can  be  defined  in  terms  of  the  directive  gain  as 


4ttU 

p 

rad 


A  very  convenient  means  of  displaying  the  results  of  the  program 
is  through  a  polar  plot  representation.  However,  because  of  the 
difficulty  of  delivering  standard  plot  routines  from  one  computer 
system  to  another,  our  plot  package  is  not  included  as  an  integral 
part  of  this  computer  code.  A  simple  polar  plot  routine  is  given 
in  Appendix  I  which  can  be  used  if  desired. 

The  next  chapter  displays  the  results  in  either  polar  or  rec¬ 
tangular  dB  plots  to  compare  against  measured  results  whenever  pos¬ 
sible.  The  results  are  normalized  to  either  0  dB  or  the  measured 
patterns  maximum. 


t 


l 


64 


thapter  VI 

APPLICATION  OF  CODE  TO  SEVERAL  EXAMPLES 

The  following  nine  examples  are  used  to  illustrate  the  various 
features  of  the  computer  code.  Each  example  is  designed  to  show 
how  a  set  of  commands  can  be  put  together  to  solve  a  single  problem 
or  a  group  of  problems.  In  most  cases,  the  input  data  sets  can  be 
constructed  in  more  than  one  way  to  accomplish  the  same  results. 

The  particular  form  of  these  examples  have  been  chosen  so  that  all 
the  commands  are  used.  As  an  aid  to  the  user,  an  echo  of  the  input 
data  is  given  on  the  line  printer  in  the  form  that  the  computer  code 
has  interpreted  the  data.  This  is  useful  for  checking  that  the  cor¬ 
rect  problem  has  been  properly  constructed.  Also,  messages  are  given 
when  the  code  misinterprets  the  data  or  when  an  error  has  been  made 
in  the  input  set.  This  makes  it  easier  to  debug  the  input  data  sets. 
Example  1A  illustrates  this  type  of  print  out.  The  other  examples 
do  not  show  this  output  in  order  to  save  space  in  this  report. 

The  computer  code  has  a  default  list  at  the  beginning  of  the 
main  program.  This  list  can  be  set  up  at  the  convenience  of  the 
irer.  If  the  defaults  are  set  up  correctly  for  the  particular  ap¬ 
plications  of  the  user,  the  same  data  will  not  have  to  be  input  in 
the  data  sets  every  time.  For  example,  the  default  list  is  set  up 
initially  to  have  the  code  give  a  line  printer  output  of  the  results. 
Since  most  user's  will  want  this  output  all  the  time  the  LP:  command 
need  never  be  specified  as  shown  in  the  examples  that  follow.  How¬ 
ever,  the  LP:  command  can  be  used  to  suppress  this  output  if  desired. 
The  pen  plotter  command,  PP: ,  on  the  other  hand,  has  been  set  false 
initially.  This  is  because  most  computer  facilities  have  different 
procedures  for  plotter  output.  Once  the  user  determines  the  best 
way  to  use  this  command  for  his  needs  and  the  appropriate  plot  routines 
are  included  in  the  code,  the  PP:  command  can  be  called  to  instate 
the  plots  or  the  default  list  can  be  changed  accordingly. 
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In  the  examples  that  follow,  all  the  results  have  been  shown 
graphically  in  some  form.  This  is  the  most  concise  way  to  show  the 
results  and  to  illustrate  the  validity  of  the  codes  operation  by 
comparing  against  measurements.  A  few  of  the  examples,  however, 
contain  the  line  printer  output  of  the  results  so  that  the  output 
numbers  can  be  checked  to  verify  that  the  computer  being  used  is 
giving  correct  results.  Different  computers  have  different  accu¬ 
racies  so  that  the  numbers  may  not  check  in  the  last  few  decimal 
places.  Example  9  can  be  conveniently  used  to  check  the  operation 
of  the  code  on  a  new  computer.  It  contains  three  examples  that  have 
line  printer  output. 


Example  1A.  Consider  the  pattern  of  an  electric  dipole  in  the 
presence  of  a  finite  ground  plane  as  shown  in  Figure  10a.  The  input 
data  for  the  H-plane  pattern  is  given  by 
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Figure  ^a.  Dipole  in  presence  of  a  square  ground  plane. 


Figure  10b.  Dipole  in  presence  of  a  square  ground  plane  with  a  pattern 
cut  corresponding  to  Example  lc. 


The  computer  code  prints  out  on  the  line  printer  the  information 
in  Figure  11  pertaining  to  the  input  data.  This  information  can 
help  the  user  decipher  how  the  computer  code  interpreted  the  input 
data.  It  also  provides  messages  to  the  user  if  the  input  data  is 
found  to  be  incorrect  by  the  code. 

The  calculated  results  for  the  electric  field  are  printed  out 
on  the  line  printer  as  shown  in  Figure  12.  The  output  shown  in  Fig¬ 
ure  12  is  for  10°  increments  from  0°  to  180°.  This  is  from  the  input 
data  in  Example  9.  Note  that  the  normalized  data  columns  may  be 
slightly  different  for  the  10°  increment  case  than  for  the  1°  in¬ 
crement  case,  since  the  maximum  value  may  be  found  to  be  a  different 
number. 

The  Eqp  pattern  is  compared  with  the  measured  results  in  Figure 
^a.  The  E^p  pattern  is  not  plotted  because  it  is  of  negligible 
value. 
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Figure  11.  Line  printer  output  for  the  code's  interpretation  of 
the  input  data  set  of  Example  1A.  The  figure  is 
continued  on  the  next  page. 
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Figure  12a.  line  printer  output  for  the  Efln  fields  of  Example  1A  and  Example  9. 
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Figure  12b.  Line  printer  output  for  the  E  fields  of  Example  1A  and  Example  9 
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Figure  12c.  Line  printer  output  for  the  total  radiation  intensity  of  Example  1A  and  Exampl 


Example  II?.  Consider  the  C-plane  pattern  of  the  electric  dipole 
in  the  presence  of  the  finite  ground  plane  in  Figure  10a.  This  prob¬ 
lem  is  the  same  as  Example  1A  except  that  the  pattern  cut  information 
is  changed  so  that  the  phi  angle  is  fixed  and  the  theta  angle  is 
varied.  The  input  data  is  given  by 


Cr*  PLAIT-  it-bi,  LXAtM'l.h  I  b. 

••  •  M  i'O  IN  I  NCI  L  • 

rid  r i:NCY  IN  UI,C. 

O  .  1 r 

I>ui  i- AiliiHIl  Cl/i 

•  i  r,  #  Vtl.  i !  • 

r  ,C . 

v‘ ,  j  ON ,  I 

PCI  PI. A  i  i:  Cl:O  b:TliY 

d 

<j  •  f  .  L «  —  •  0 

0  •  f  "s.j  .  b  9  b  •  b 

.  b  « “  o  •  b 
■) .  i  j  *  t  f  ™  J  . 

but  bOUivCi-  (JirCdiuPY 

-)  •  1  P  9  Ir  •  9  IT  . 

l)  .  I  I'.  |  V  I  *  .  f  V  • 

N  9  C  •  1)  9  l  9 

I 

.*■  ‘  LXPbVil;  CO'. i|; 

tl  :  CLOP 

The  E0p  pattern  is  compared  with  the  measured  results  in  Figure 
13b.  The  E^  pattern  is  not  plotted  because  it  is  of  negligible 
value. 
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MEASURED 

CALCULATED 


Example  iC.  Consider  the  pattern  of  the  electric  dipole  taken 
across  the  corner  of  the  finite  ground  plane  in  Fioure  10b.  This 
problem  is  the  same  as  in  Example  1A  except  that  the  pattern  cut 
coordinate  system  is  changed.  The  input  data  is  given  by 


l  *  :  l.Aii  t «  i  -  X  A  I  >.  . 

*  • : r-i i  t ; >  i r:  I  'cid  s 
_> 

1-tC  r  HcOl.i.!  it:  Y  IN  Oil/., 
o .  t 

1  A;  it-.iJh  CiJ  i 

•Oj.  .  ,yj.  ,o. 

i ,  s  v . . 

0  ,  ~ ,  I 

.-o:  1-LA  1 1:  UbUihTKY 

4 

■ !  .  .  —  j  *  )  |  J 
4  .  (  .  )  |  “  J  . 

’*•  *••.  b  y  —  J.  • 1  > 

LUur<Ci:  (ji:0;.’,l:TT<Y 

*).  1  ^  i  i  .  f  • 

■4  .  y  " 1  .  (S0.  y  l  . 

*  r  .l)  il-  . 

*  .  ^  • 

OS  LXcCP'if:  C:0|)l: 

'  (J  <Jl  l>|- 


The  Ea  pattern  is  compared  with  its  measured  result  in  Figure 

13c  and  the  F.  pattern  compared  with  its  measured  result  in  Figure 
<pP 

13d. 
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Example  2.  Consider  the  pattern  of  a  A/2  slot  antenna  mounted 
in  the  center  of  a  square  ground  plane  as  shown  in  Figure  14.  The 
pattern  cut  is  taken  across  the  corner  of  the  plate.  This  example 
illustrates  how  the  TO:  command  can  be  used  to  show  the  strength 
of  a  particular  diffraction  mechanism  by  removing  the  mechanism  from 
the  computation  for  comparison.  In  this  case  the  corner  diffracted 
fields  are  removed  in  the  second  execution  by  setting  LCORNR*. FALSE. 
Of  course,  other  fields  can  be  modified  in  other  problems  by  using 
the  JMN's  and  JMX's  numbers.  The  input  data  for  this  example  is 
given  by 


a-.*  PLATE  iEST,  EXAMPLE  2. 

UN  * 

J 

i-KS 

I  0. 

PU*  PATTERN  CUT  OVER  CORNER 

0  .  ,  O  •  ,  S  0  .  ,  it  *  . 

E,4  0. 

0 , COO',  I 

PCs 

4 

6 . ,  6.  , 0 . 

-o.  |  — c. , o  • 

o. , “O. ,0. 

CO* 

0. , 0.  ,0  . 

90.  ,V0.  ,90.  ,0. 

I ,0.0,0. 

I  .  ,0. 

XU* 

10*  REMOVE  CORNER  OIEHRACTI ON 
E ,  r  ,  r 
i  ,k  ,r- 

1,7,1  ,3,1,4 
XU* 

EN* 

The  results  with  and  without  the  corner  diffracted  fields  in¬ 
cluded  are  compared  in  Figure  15a  and  b  for  the  and  E.„  fields, 

ep  4>p  ’ 

respectively. 


-  WITH  CORNER 

DIFFRACTION 

--  WITHOUT  CORNER 
DIFFRACTION 


? 

> 

V 


Figure  15b.  Comparison  of  E.  pattern  with  and  without  corner 
diffracted  fiela5  for  a  half-wave  slot  antenna 
mounted  on  a  square  plate. 
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Example  3.  Consider  the  pattern  of  an  electric  dipole  in  the 
presence  of  an  eight-sided  box  as  shown  in  Figure  16.  The  input  data 
is  given  by 


Cfc*  blUHi  SIDED  10X  TEST,  EXAMPLE  3. 

t-Ht 

V.y4 

PL)* 

id  . ,  0 .  ,  V  V)  .  ,  0  . 

1,V0. 
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SG* 
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—  .  1  22  ,  —  .  1  02  3 , .  1 
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— .  122,.  1 023  ,  - .  I 

PG*  U:hi  HACK 
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0. ,-. 1707 ,-. 1 
0.,-.  1707, .  1 

—  .122,—  .1023,.! 

-.1 22  ,-.1023,-.  I 
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The  E„„  pattern  is  compared  with  measured  results  fill  in  Figure 
The  £  pattern  is  not  plotted  because  it  is  of  negligible  value. 
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electric  dipole  in  the  presence  of  an  eight  sided  box. 


Example  4.  Consider  an  electric  dipole  in  the  presence  of  a 
finite  circular  cylinder  as  shown  in  Figure  18.  This  example  illus¬ 
trates  how  once  the  cylinder  geometry  is  set  up  different  cases  can 
be  run  by  just  varying  the  data  which  needs  to  be  changed  for  that 
particular  case.  The  input  data  is  given  by 


Cl- 1  CYLINDER  TES) .  EXAMPLE  4A. 

fi< » 

y.y* 

pu* 

0  .  ,  0  •  f  VW.  ,  0  . 
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vo. ,w.t 180. ,0. 
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PDi  CHANGE  PATTERN  CUT 

0. , 0. ,V0. ,0. 
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XC* 

CL*  CYLINDER  TEST,  EXAMPLE  4C. 

NS*  CALL  FOR  NEW  SOURCE 

SG» 

0.0"/o,0.  ,0.2 

V0. ,0., 180. ,0. 

0,0.8  ,0. 

1  .  ,0. 

aU* 

;.*?•  * 


The  line  printer  output  of  the  results  for  Example  4A  are  shown 
in  Figure  19.  The  input  data  is  really  that  of  Example  9  where  only 
half  the  pattern  is  computed  in  10°  increments. 
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Fi^jrp  19a.  Unp  printer  output  for  the  E  fields  of  Example  4  and  Example  9 


AO-A097  416  OHIO  STATE  UNIV  COLUMBUS  ELECTROSCIENCE  LAB  F/G  20/14 

NUMERICAL  ELECTROMAGNETIC  CODE  (NEC)-BASIC  SCATTERING  CODE.  PAR— ETOUl 
SEP  79  R  J  MARHEFKA.  W  D  BURNSIOE  N00123-76-C-1371 

UNCLASSIFIED  ESL-784508-18  NL 


Figure  19b.  Line  printer  output  for  the  E.  fields  of  Example  4A  and  Example  9 


gurp  1Qc.  Line  printer  output,  for  the  total  radiation  intensity  of  Example  4A  and  Example 


The  calculated  results  are  compared  with  measured  resultsflll 
in  the  following  figures. 

The  E^  pattern  for  Example  4A  is  shown  in  Figure  20a. 

The  E^p  pattern  for  Example  4B  is  shown  in  Figure  20b. 

The  E^  pattern  for  Example  4C  is  shown  in  Figure  20c. 

The  EeP  pattern  for  all  three  cases  are  not  shown  because  they 

are  of  negligible  value. 
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der  (pattern 


measured  (Bach)  and  calculated  radiation  pattern  for  ExD 
the  x-z  plane  parallel  to  the  x-axis  in  the  presence  or 


Example  5.  Consider  the  pattern  of  a  magnetic  dipole  in  the 
presence  of  an  elliptic  cylinder  as  shown  in  Figure  21.  This  example 
illustrates  how  the  LOUT  parameter  in  the  TO:  command  can  be  used 
to  print  out  the  individual  fields  reflected  and  diffracted  by  the 
body  under  consideration.  Note,  also,  that  since  the  units  and  fre¬ 
quency  are  not  specified  in  the  input  set  the  input  is  therefore 
assumed  to  be  given  in  wavelengths.  The  input  data  is  given  by 


CE»  tiLJ.il  TiC  CYLINDER  TEST,  EXAMPLE  b . 

Pi  >  * 

0  .  ,  •  ,  90 •  ,  i>;  . 

i  ,90. 
v> ,  3  00 ,  I 
cut 

0  .  ,  I  . 

-000 .  ,  90.  ,  ECU  .  ,  yn. 

SO  • 

l  .000  ,0  .»0L,O. 

i 0  •  ,  0  •  ,90. ,0. 

I  ,  hi  .  J  , '/)  . 

I  .  ,  0. 

XCt 

10*  PR  INI  I  HI '  I  V  I  DUAL  FIELDS  AROUND  EHADOR  BOUNDARY 
F , r  ,  i 
i  .  i  . I- 

Hi;*  REDUCE  ANGULAR  RANGE 

0. ,0. ,90. ,u. 

1,90. 

0O<  ,009  ,  I 
AO* 

E  t,  * 


The  reflected  and  diffracted  fields  in  the  region  close  to  one 
of  the  shadow  boundaries  of  the  elliptic  cylinder,  as  printed  by 
the  line  printer,  are  shown  in  Figure  22.  The  different  types  of 
fields  can  be  interpreted  by  looking  up  the  integer  indices  in  Table 
3.  The  first  two  columns  of  real  numbers  are  the  magnitude  and  phase 
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of  the  field  and  the  second  two  columns  of  real  numbers  are  the 
magnitude  and  phase  of  the  field.  The  polarization  is  referred 
to  the  reference  coordinate  system  for  this  type  of  print  out. 

The  E^p  pattern  is  plotted  in  Figure  23  compared  with  a  moment 
method  solution.  The  EQp  pattern  is  not  plotted  because  it  is  of  neg¬ 
ligible  value. 


Figure  21.  Elliptic  cylinder  configuration  excited  by  a 
magnetic  source  parallel  to  the  z-axis. 


cnattii  v.’rtT-  rrr/t‘-‘-  ooconr 


s- 

<L> 

■o 


im  tt>  < x  8«e  »  »  *  r  <  tf  ® 

.*  fr^oo4^v~vf'^^«»voc»oooo<rr-«qB**^»'0^^iA'^^ 
>oao-<ofM<-<^v'»«‘^r^oo-^^€aooooo«'«NO'^H»'CN-^-^®r-o^ 

9  9  90000000P»HN*«^-<Nf9a 
rt'n(M9^mi«fl«09  <9  O  O  C*  9-  OOOOOOlT^©r~  X<  '•■>  O  fy  **  O  oa 


r  ry  /n  &>  tr  ^-*<Nvfl^M«n»»oo' 

it  »  | 

-*  tit  •— *  •  t  *  I  1  I  I 

1  I  I 


o  tr>  ®  r»  r-  r»  |  _<  ry  m  ^ 
I  I  I  I  I  I  I  (  t  I  I 


«»'.'NoiraDoomiroomo*-4'^P,“r‘-4^c>ooooooooooo0^o^«f^^ 

9<p>'trir-<rac>oir^^><vvcsc9'7i/‘irooooooooooootro«r^9'ir 

.— 4  »T  «T  — 4  ,£-  <y  ~4  r-«  «•  H  N  N  «  000000000000'9',»'>«^«N^ 

IT  3  V«a9iTirv^«vp<£  C^XOOOOOOOOOOOOiTSi^^^^ 

-«t/*oir  CMroif  o»r^-«-<r-«^«-4-4ooooooooOOoooOi-«^^^ 

OUM-JOOHHOO^H^HHH-JOOO.^OOOSOOOOHHH^HM 


r'*^'9or-«>*<Sa?r»-^-^9-^oooDoeOooooO*9^CO'^^'*'^®r^p,“<^^ 
lro-^-^^r-f^'^f'ir^<r-r'^^c»ooooooo^^'*>ir'r^O‘^Pi(ncrr,<W' 
'^NfNof*^9ir(ir(r(AO''('i«ooooooooir^^®N*Hir^,>,vNp4 
NiTif'nwo9B5^f^"*9r  r*r>^^c»ooooo^f**»r9''^®'«,^»r'»',,^« 
iA®«5jr^99^r‘mp3«**ic«<*'"OOOOOON^CB®'«or<»®««o 
a'W3'«T»fi'ni/‘«rtrn  —  x  ^-oooooo^f'-r»»^«r-y«9f'-f^'^^'y 

^»l  |r-f  ««« 

-♦  -•  I  t  I  I  »  ^  ^ 

I  I  t  I  I  I  t  t  t  |  f 


OoOOOoOOOOOO^OOOoOOOOOOOOOC>OOOOoOOO 

oooooooooooor'oooooooooooootsoooooooo 

©000©000©000©0000000000000©0©©00000 

ooooe»©oooooo.-»oo©oooooooooo©o©ooo©©© 

0000©000000000000000000000900©00000 
■  0000900o0000^0000000000000<900000000 


000-4000-4000-40 


I  ♦  IT  ^  N  •  9 

o  ©  o  ©  ©  © 

ry  ry  ry  ry  ry  ry 


looyoooi/'  ooo  x>  o  "■  ooc  o  9  y  vr  r  n  id  9  *  /  f  r*  ®  9  v  it  ^  ®  9 

©  O  o  7*  O  OTOOt>OOOOOOOOC>OPOOO 

ry  ry  ry  -;  ry  ry  ry  ry  ry  ry  ry  . .:  «ryryryryryryryryryryry 

nooyo^oirnootO^OBOtyrf1  ^  r-  <*>  9.  ^  tr  «r>09oooooo 
yy»*io-4N*io-«»y',l'o"'">*io^o30ooooooooooe»oo{>oo 
y^HN-«^-<'y-»^^<,y^N-4n*-«(Ny(Ny'NnM<yfi(Ni,yiyNOOoooo 


>> 

U  t. 

a> 

<J  -O 
•*“  E 

+->  3 
a  c  *o 

•r-  C 
r-  X  «0 

r-  a; 

0)13  0) 

C  T3 
Q)  •«-  3 

_c  +j 

P  S-*r 

aj  c 
>>  cn  a> 

-o  a>  <0 
_  +-»  E 
T3  c 

<u  -r-  a> 

«M  _c: 

U  Ol  «M 
*3  -C 

i-  +->  CO 
<4-  +-> 

^  c 
•f-  o  a) 

"O  t/)  • 

C  (U  >i 
to  O  V-  «— 

-a  -r-  Ql  cu 
r—  4->  <d  > 
<U  *r~  l-  ‘r- 
•r-  C  -M 
4-  **■“  tO  V 
M—  i-  <D 

»—  a>  a;  o. 
<tj  "o  jo  l/> 
3  E  a> 
■o  a>  3  s- 
•r-  -c  c 
>  h- 

•*-  • —  t/> 

■O  fO  TD 

c  •  ai  *— 

•|“  UJ  1-  <u 
ZD  •«“ 

cu  o:  <u  d- 


11  -©■ 

4-  h-  LU 

O  Z3  . 

O  CO  "O 
-4->  — I  C 

3  0)  #0 

Q--C  r— 

•M  «•->  XI  Q. 
3  •»-  *3  O 
Oil —  LU 

Lin  c  a> 

0)  *r*  _c 

-m  a;  -m 

c  r-  c 
•*-  o.  a>  <♦- 
v.  E  >  o 
Q.  «J  •»— 
x  cnai 

QJ  LU  tO 

c  a>  03 
4-  .rr 
_j  o  hi  a 


100 


- GTD 

•  •  •  MOMENT 
METHOD 


Example  6.  Consider  an  electric  dipole  in  the  presence  of  a 

plate  and  a  finite  circular  cylinder  as  shown  in  Figure  24.  This 

example  illustrates  how  the  input  data  can  be  manipulated  to  analyze 

the  effects  of  scattering  bodies  separately  and  in  combination  with 

one  another.  The  following  input  is  shown  as  if  all  four  cases  are 

run  consecutively  in  one  run.  Of  course,  the  input  could  easily 

be  constructed  for  individual  runs  for  each  case. 

Chi  PLATE-  AND  CYLINDER  TEST,  EXAMPLE  6. 

PD  i 

0 .  ,  k).  ,VW.  ,  w . 

i  ,  Vt). 

U  ,  ..  Ok'1  ,  I 
UN  i 

j 

Nil 

4  . 

DC  i 

vJ .  ,  b .  CTb,W. 

Vv). ,  IP  • , 0 • , O • 

W  ,  IrJ  .  "V  ,  t)  .  Ip 

I  .  «  tv . 

PCI 

4 

b .  ,  IP.  ,‘j  . 
b  •  ,  t  •  •  ,  —  tv . 

-b.  ,IP  .  ,  -b  . 

-b.  « k > .  , b. 

rt  ’l  * 

v).  ,  IP  .o  12b, -0.  12b 
o.  ,c.  ,ct^.  ,i. . 

CD  * 

I  .2b,  1 .2m 

-E.  b,  VC- .  ,  ri.  b  ,  VD . 

ACM  b OUNCE,  PLATE  AND  CYLINDER  TEST 

DPI  REMOVE  PLATE 

XO*  SOUkCE  AND  CYLINDER  TEST 

NCI  R I-j.IOV h  CYLINDER 

XO I  SOURCE  IE Si 

Hi  I  Rl  lNLi  1I.I/.E  ORIGIN 

0 • ,  t > •  ,0. 
ip  •  ,  i  /  •  ,  V  O  •  ,  t  . 

Pi  i  i  it  ED  hr  I  N  E  PLAYE 

4 

ib. 

,n.  ,-b. 

~  *  ,  t ;  .  ,  — b  . 

*”  J.  ,  P  • ,  b . 

m.s  SOURCI:  AND  PLATE  TEST 
i  r.s 


102 


The  line  printer  output  of  the  results  for  the  source,  plate 
and  cylinder  in  combination  are  shown  in  Figure  25.  The  input  data 
is  really  that  of  Example  9  where  only  half  the  pattern  is  computed 
in  10°  increments. 

The  calculated  E^  patterns  are  compared  with  measured  results 
in  the  following  figures. 

The  E^  pattern  for  the  source  alone  is  shown  in  Figure  26a. 
The  E^p  pattern  for  the  source  and  plate  is  shown  in  Figure 

26b. 

The  E^p  pattern  for  the  source  and  cylinder  is  shown  in  Figure 

26c. 

The  E^p  pattern  for  the  source,  plate  and  cylinder  is  shown 
in  Figure  26d. 

The  E0p  pattern  for  all  four  cases  are  not  shown  because  they 
are  of  negligible  value. 
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Figure  25b.  Line  printer  output  for  the  E  fields  of  Example  6  and  Example  9. 


:igure  25c.  Line  printer  output  for  the  total  radiation  intensity  of  Example  6  and  Example 
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Example  7.  Consider  a  slot  mounted  on  the  wing  of  a  Boeing 

737  aircraft.  The  computer  model  of  the  Boeing  737  is  illustrated 

in  Figure  27.  The  input  data  is  given  by 

vj.'.'i  /  Im.'i.Ari  i-XntPLF  7. 
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The  E^p  pattern  is  compared  with  its  measured  result  in  Figure 
28.  The  measurement  was  made  on  a  1/20  scale  model  of  a  Boeing  737 
at  NASA  (Hampton,  Virgina).  The  antenna  is  a  KA  band  waveguide  mounted 
in  the  wing[8].  The  Eq  pattern  is  significant  for  this  case,  how¬ 
ever,  it  is  not  shown. 


Comparison  of  measured  and  calculated  E„  results. 
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Example  8A.  Consider  the  pattern  of  a  set  of  four  electric 
dipoles  in  the  presence  of  a  square  plate  over  an  infinite  ground 
plane  as  shown  in  Figure  ?9.  The  currents  are  specified  by  the  NEC- 
Moment  Method  Code  as  used  at  NOSC.  The  input  data  is  given  by 
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f  £_  J  wOfS 

.  3bbC0 

.20637 

.29404 

. JbbOO 

.20637 

. 2bboO 

. 3bbOO 

.20637 

. 4 1 ob6 

. 3b56G 

.20637 

. 4 7 7b2 

. 3bb6fcl 

.20637 

.32bbl:-3 

. 1 796b— 2 

•Obc4H-3 

. 4092 E-2 

•  7  9  7 3H-3 

.  bo2  i  H-2 

.6  7  obi: -2 

.  40H9E-2 

.3462H-2 

.  r/y 4f?-2 

•  3  4o2h-3 

.  1794H-2 

.o73bi:-2 

.  4089 H-2 

. 7  C7 3H-2 

. b62 1 H-2 

.obo4t:-3 

.4092 H-2 

. o2bbH-2 

.  17 90 H-2 

.32bbt-2 

.  1 7 90 H-2 

.obo4b-3 

.4 092 H-2 

.  7  C7  31: -2 

.  b62 1 H-2 

•o73bh-3 

.4 OHS H-2 

.  34o2H— o 

.  17 94 H-2 

.04O2H-2 

.  1  7o<  F-2 

. o  7 j  bh-3 

.4089 H-2 

.  / C/7 3H-3 

. bo2 1 H-2 

•  O  bo  4b— 3 

.  4092  H—  2 

.o2H8b-3 

.  I7VOH-2 

X()« 

NG«  KbMOVL- 

GHOtlM)  PLAN 

H 

XU  * 
nN « 


. 06096 

.  00000 

. 00000 

. 06096 

.  00  000 

.  00000 

. 06O96 

.  00000 

. 00000 

. 06096 

.  00000 

.011(100 

.06096 

.  00(100 

. 00000 

. 06096 

. 00000 

. 00000 

. 0009 6 

. 00000 

. 00000 

.06096 

. ooooo 

. 00000 

. 06090 

.  00000 

. 00000 

. 06090 

. 00000 

. 00000 

. 06096 

.00000 

. 00000 

.060196 

.  0000(1 

. 0000" 

.06090 

. 00000 

. 00000 

.06096 

. 00000 

. 00000 

.06090 

.  00,000 

. 00000 

.06090 

. 00000 

. 00000 

.06096 

. 000/0 

. 00000 

.06096 

. 00000 

. 00000 

.06096 

. 0O000 

. 00000 

. 06096 

. 00000 

. 00000 
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The  line  printer  output  for  the  first  execution  with  the  in¬ 
finite  ground  plane  is  shown  in  Figure  30.  The  results  are  shown 
for  10°  increments  for  brevity.  Since  a  range  is  specified  the  fields 
are  in  volts/meter.  Note  that  the  range  is  in  the  far  field  of  the 

p 

maximum  dimension  of  the  plate  (R>2D  /A).  Also,  since  the  power 
radiated  is  known  from  the  NEC-Moment  Method  Code  the  power  results 
are  given  in  terms  of  directive  gain  rather  than  in  terms  of  the 
radiation  intensity. 

The  directive  gain  normalized  to  isotropic  is  plotted  in  Figure 
31a  for  the  <t>=0°  plane.  The  result  is  compared  against  the  infinite 
ground  plane  case  and  the  case  for  the  plate  in  free  space  which 
is  given  as  the  second  execution  in  the  input  set.  Similarly,  the 
directive  gains  for  the  three  different  cases  are  plotted  in  Figure 
31b  for  the  4>=90°  plane. 


Figure  ?9.  Geometry  for  the  problems  of  dipoles  over  a  square  plate  and 
infinite  ground  plane  showing  the  side  and  top  view. 
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Figure  31b.  Comparison  of  the  directive  gain  of  four  dipoles  over 
an  infinite  ground  with  four  dipoles  over  a  square 
plate  over  an  infinite  ground  and  four  dipoles 
over  a  square  plate  alone  (4>=90  , 
horizontal  polarization). 
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Example  8B.  Consider  the  pattern  of  a  set  of  four  electric 


dipoles  in  the  presence  of  a  square  plate  over 
plane  as  shown  in  Figure  29.  The  four  dipoles 
analytic  representation  using  the  SG:  command, 
the  use  of  the  PR:  and  US:  command.  The  input 


an  infinite  ground 
are  specified  by  their 
This  example  illustrates 
data  is  given  by 


Cb*  N bC  ibSI,  bXAMPLL  8H. 
PD» 

0  .  ,  0  .  ,  V  Id  .  ,  0  . 

b 

0,90,  I 
bK» 

0 .29vb 
UM 

3 

HG» 

1  bU . 

UP  J 
K  1  s 

0.  ,  H.  ,2b. b 
0. , 0. ,90. ,0. 

PGs 

4 

2b. ,-2b . ,0. 

2b. ,2b. ,0. 

~2b  .  ,2b  .  ,0. 

-2b . , -2b. ,0. 

UM 

1 

USs 

I 

PM 

. 00  I j3b 

SG«  SCURCb  #1 
- .obbc0  ,-. 3bb60 ,  .20637 

90. . 0.. 90.,  90. 

0,0.3048,0. 

. 00bo73 ,82.2 
SG*  SOURGb  #2 
. 3bbo0 ,-. 3'jbot'',  .20037 

90. . 0.. 90. .90. 

0,0.3048,0. 

.00 bo ',3 ,82. 2 
SGt  SOUWCb  #3 

- .3bjo0, . 3bbo0, .20637 

90. ,0., 90. ,90. 

0,u. 3048,0. 

. 00bo7J ,82.2 
SGI  SOJHC1:  #4 
. 3bbO0, .3bb60,. 20627 

90. ,0., 90., 90. 

0,0.3048,0. 

. 00 bo 7  3 ,82.2 
XU  * 
i:N  I 
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This  example  gives  approximately  the  same  results  as  those  in 
Example  8A.  The  directive  gain  is  within  approximately  0.5  dB  through¬ 
out  the  entire  pattern  range.  The  pattern  shape  is  essentially  identi¬ 
cal  to  the  one  in  Figure  30a  for  the  four  dipoles  over  a  plate  over  an 
infinite  ground  plane  so  the  pattern  is  not  shown. 


1  ?5 


Example  9.  Consider  a  combination  of  three  different  problems. 
The  first  problem  is  the  plate  test  of  Example  1A.  The  second  problem 
is  the  cylinder  test  of  Example  4A,  and  the  third  example  is  the 
plate  and  cylinder  test  of  Example  6.  This  example  illustrates  the 
use  of  the  NX:  command  in  defining  entirely  different  problems  in 
the  same  input  set.  This  is  a  handy  example  that  can  be  used  for 
the  initial  start  up  test  for  the  code  on  a  new  computer  system, 
since  the  plate  and  cylinder  are  tested  separately  and  in  combination 
with  one  another  in  the  same  input  set.  The  input  data  is  given 
by 

CM*  MULTIPLE  TEST  EXAMPLE 

Cfc  *  PLATE  iESi,  EXAMPLE  I  A. 

'JN* 

J 

Elf  * 

Pc* 

0  •  ,  kv  •  ,  V  0  .  ,  w  . 

0,  iHki,  in 

PCX 

4 

0  • ,  ***3  •  b  ,  3  •  b 
VI .  ,  -3  .b  ,-3.  5 
0  . ,  3  .  S,  -3  .  b 

sc* 

b  •  I  0  ,  V  •  ,  0  . 

0  •  ,  xi  m  ,  V0  .  ,  0  • 

0 , .  b ,  v . 

I  . ,  c  • 

AO  * 

NX* 

CE*  CYLINDER  TEST,  EXAMPLE  4A. 

EP* 

V.  V4 
Pl>« 

0.  ,  lr}.  ,V0.  ,0  . 

1  »V0. 

*0,3  /0,  10 
SC* 

0  .  ,  0  .  I  S*  ,  0  • 

*0.  ,0.  ,  lfcJ0.  ,0. 

0 , 0  •  ,  0  » 

I  .  ,  V  *  . 

cc* 

0  .  I  ,  V)  •  I 
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-fcl.  I  I  ,VM.  ,10.  I  I  ,  VM. 

XU* 

NX* 

i;i:i  PLA r*  AND  CVI.l NDhli  TEST,  EXAMPLE  6. 
PI;  s 

D  .  ,D.  ,S>D.  ,l'  . 

1  ,  S4). 

VM,2  /t,,  K) 

ON* 

.0 

r  l<  * 

4  . 

SO* 

vO.  ,  b.  O2b,l0. 

vtf.  ,u  . ,  to.  ,fcO . 

vMO.bLb.H.V; 

I  .  ,H. 

P(j* 

4 

b  .  ,  14 .  ,  b  . 
b . ,  v) .  ,  -  b . 

~b .  ,  id  . ,  ™  b . 

-b.  ,») . ,  b. 

K  1  * 

>4.  ,  12  12b, -(4.  12b 

.  ,2.  ,  V0 .  ,14  . 

CU  * 

I  .2b,  I .2b 

-H.b,U0.,d.b,VM. 

XU* 

Eh* 


The  line  printer  output  and  the  graphical  representation  of 
the  fields  for  each  of  the  three  examples  above  are  given  in  their 
respective  example  sections.  The  line  printer  output  of  the  first 
execution  is  given  in  Figure  12  and  the  plotted  output  is  given  in 
Figure  13a.  The  line  printer  output  of  the  second  execution  is  given 
in  Figure  19  and  the  plotted  output  is  given  in  Figure  20a.  The 
line  printer  output  of  the  last  execution  is  given  in  Figure  25  and 
the  plotted  output  is  given  in  Figure  26d. 
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The  following  is  a  listing  of  a  polar  plot  subroutine  which 
can  be  used  to  generate  the  polar  patterns  given  in  Chapter  VI.  Note 
that  this  code  refers  to  two  subroutines  "PLOTS"  and  "PLOT"  which 
must  be  added  by  the  system  or  operator.  The  definitions  of  these 
routines  are  given  in  the  comments  associated  with  the  code. 


SUBROUTINE  POLDLT ( ET,  RP.  I  PLT,  IPHS.IDM) 

T‘>IS  1.  OUT  IMF  IS  USED  jo  Pf.nl  THH  RESULTS  If  TERMS  OF  A 
• ' n I . ft R  PLOT.  TllF  CAM.  TO  SUBROUTINE  "PLOTS"  IS  "SED  TO  INITIALIZE 
i  ■*:;  M.OlTEN.  T1JF  CALLS  TO  TME  SUBROUTINE  "PLOT"  ARE  USED  TO 
DRAW  ■*»■:?-  LATA  AS  FOLLOWS* 

CALL  PLOT ( X , Y, N) 

X,Y=COOFm NATES  OF  THE  NEW  PLOT  POINT. 

N=2  P- T'  15  OOWN  MOV  I  MO  TO  THE  N  F’N  POINT. 

L=3  PEN  IS  »JP  MOVING  TO  TMF.  *»EW  POINT. 

N=9W  F-UFFER  USED  TO  STORE  PLOT  DATA  IS  EMPTIED  TO  PLOTTER 

N«0  IMPLIES  ORIGIN  SHIFT  TO  THE  *’EW  Pf)IMT. 

N>0  IMPLIES  NO  ORIGIN  SHIFT  AFTER  MOVING  TO  NEW  POINT. 

COMPLEX  HT< I  DM) 

DIMENSION  IBUF(ICC) 

DATA  ni  ,TPI  ,l)PR/3.  I  a  I  59265, 6.  2831853 , 57  .?9b71'?b3/ 
r.  MX-O. 

DO  101  IP=0,560, I PHS 
I=IP+I 

f:M=CABS(FT(  I  )) 

I  F(  EfA.GT .  EMX  )  EMX=E» 

CONTINUE 

CALL  uLOTS( I BUF, 100,3) 

CALL  PLOT (4.25,5.S,-3  ) 

POLAR  GRID  **♦ 

DO  110  1  =  1,^ 

RG=RP*I /4 . 

CALL  PL0T(RG,O.  ,3  ) 

H)  I  1  O  J=D. 360. 2 
ANG=J/i)PR 
>.  X=W(  *COS  (  AN6  ) 

YY=R(  *SIN< AMG) 

CALL  PLOT (XX.YY.2) 

DO  112  1=1,6 
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a r.’u=(  i- 1  )*pi/6. 

ANGS=ANO-*-pI 

ANGF=ANG 

I  F<  I  .  FO.  ?*(  J/P  )  )  GO  ;•  t  1 
A  NOS=  AN’G 
A  NGr  =  AfJG+PI 
111  CONTINUE 

AX=RP*COS(AM(,5  ) 

V  Y=HP*SI N  (  ANGO ) 

CALL  PLOT! > X , VY  .3 ) 
XX=RP*COS( ANCF) 

YY=RP*S I N ( ANCr  ) 

II?  CALL  PLOf(XX,YY,2  ) 

*  **  PATTrl'N  ^ I .: >  1  *** 

PO  I2C  IP*f).3fsC,  I  PHS 
I  =  f  P+  I 

H1M-CAP^(LT(  I  ) ) /c MX 
I FC I PLT-2  )  121, 122, 12  3 
12  1  <?r>=Rp*?-TM 

CO  TO  I 25 

122  RP»HF*t'TM*STM 
GO  TO  125 

12  3  I  F  (  FI M .  L  i'  .0.01)  FTM=0 . 0  1 
h  D=2  C . *  A  LOG  I C ( r TM ) 

I F( FO.LT. -^0. )  R0=-40. 
RP=RP*(i<n+40.  )/^0. 

125  CONTINUE 

AM0=I  P/f'PR 

X  X =.-ii  ’  *S  I M  ( AI.'C- ) 

YY=R!)*COS(ANG) 

I P1£N  =  2 

I F( I .HO. I )  IPbN=3 
1 2 0  CALL  PLOT(XX, YY, I  PEN) 

CALL  PLOT (4.25 , -5 .5,-3) 

130  CONTINUE 

CALL  pL0T(0. , 0. , ) 

RETURN' 

END 
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